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SECTION  1 
INTRODUCTION 


A problem  of  continual  concern  to  the  Air  Force  relates  to  determining  how  new 
and  existing  aircraft  affect  airfield  pavements.  Two  questions  are  usually 
asked:  (1)  can  a pavement  support  the  aircraft?  and  (2)  what  is  the  remaining 
service  life  of  the  pavement?  The  former  question  concerns  the  assigned  mis- 
sion capability  of  a given  facility,  that  is,  can  a given  pavement  safely  sup- 
port an  assigned  aircraft,  the  B- 52  for  example?  The  latter  question  concerns 
the  number  of  times  the  assigned  aircraft  can  safely  operate  from  the  pavement 
in  question. 

Historically,  two  different  techniques  have  been  used  to  answer  these  questions-- 
namely,  empirical  relationships  developed  through  airfield  performance  studies, 
and  full-scale  testing  of  pavement  sections.  Both  of  these  techniques,  however, 
have  serious  deficiencies  for  a user  who  is  required  to  maintain  airfield  pave- 
ments under  a wide  variety  of  geotechnical,  geographical,  and  climatological 
conditions.  Additionally,  destructive  testing  is  inherent  in  both  of  these 
methods,  and  thus  runways  must  be  closed  for  extended  periods  of  time  while 
material  samples  are  gathered  for  laboratory  testing.  This  process,  besides 
being  time  consuming,  has  a significant  impact  on  facility  use  and  requires 
considerable  expenditures  for  testing  and  analysis;  and  at  best  the  underly- 
ing theory  is  questionable,  though  the  results  are  conservative. 

To  improve  upon  the  situation  the  Air  Force  planned  a pavement  research  pro- 
gram to  develop  a rational  pavement  evaluation  procedure  (1967-69).  Central 
to  this  procedure  was  the  use  of  nondestructive  testing  techniques  to  obtain 
the  in-situ  properties  of  the  paving  materials.  The  advantages  of  such  an  ap- 
proach are  obvious:  test  pits  do  not  have  to  be  placed  in  the  pavement,  more 

data  can  be  obtained  during  equivalent  periods  of  runway  closure  thus  providing 
a better  measure  of  pavement  material  characteristics,  and  testing  can  be  done 
at  night  so  runways  remain  open  for  use  during  the  day. 

This  report  provides  documentation  on  much  of  the  Air  Force  sponsored  research 
directed  toward  developing  a nondestructive  pavement  evaluation  procedure.  It 


also  serves  as  a user's  manual  on  the  CERF-developed  vibratory  testing  tech- 
nique for  pavement  evaluation  studies.  The  equipment  and  analytical  tech- 
niques used  in  this  procedure  have  been  transferred  (1975)  to  the  Air  Force 
Civil  Engineering  Center  for  use  in  their  pavement  evaluation  program.  Cur- 
rent research  is  concerned  with  making  the  nondestructive  evaluation  equipment 
air- transportable. 


SECTION  2 
HISTORICAL  REVIEW 


INTRODUCTION 

This  section  presents  a historical  review  of  the  development  of  the  Air  Force's 
nondestructive  pavement  test  (NDPT)  procedure.  The  need  for  this  procedure  co- 
incides with  the  introduction  of  wide-bodied,  jumbo  jet  aircraft  (1967).  The 
research  described  in  the  following  paranraphs  does  not  include  all  of  the  pave- 
ment research  sponsored  by  AFWL ; only  those  projects  which  relate  directly  to 
the  NDPT  project  are  covered. 


RESEARCH  PROGRAM  (1967-1969) 

AFWL  outlined  a study  to  evolve  a rational  pavement  evaluation  procedure  having 
following  capabilities: 

(1)  A technique  by  which  the  load-supporting  capacity  of  an  existino  air- 
field pavement  could  be  determined.  The  physical  parameters  required 
for  this  evaluation  were  to  be  determined  by  nondestructive  testing 
techniques. 

(2)  A technique  by  which  the  stresses  and  strains  at  a general  point  with- 
in the  pavement  could  be  determined  when  the  load  is  applied  through 
the  various  gear  configurations  of  modern  aircraft. 

(3)  An  analytical  procedure  by  which  the  influence  of  the  strenoth  and 
thickness  of  the  various  pavement  components  could  be  compared  in 
terms  of  the  stress  and  strain  response  of  the  pavement  system. 

(4)  An  analytical  means  of  evaluatinq  the  useful  life  of  an  airfield 
navement  subjected  to  repetitive  loading  and  mixed  aircraft  traffic. 

(5)  Applicability  to  both  rigid  and  flexible  pavements,  composite  pave- 
ments, and  pavement  overlays. 


This  study  was  approached  through  a critical  review  of  the  literature  concerned 
with  pavement  desigi  «:nd  evaluation  procedures,  a review  of  the  new  computer 
techniques  (finite  elements)  which  were  emerging  at  that  time,  and  a review  of 
the  new  approaches  to  material  modeling  (constitutive  equations)  which  were  be- 
ing used  in  allied  fields  of  science.  The  results  of  this  study  were  published 
in  a two-volume  report  (ref.  1).  Volume  I contained  a summary  of  flexible  and 
rigid  pavement  design  methods  for  highways  and  airfields  and  an  extensive  dis- 
cussion related  to  the  theoretical  approaches  to  pavement  evaluation,  in- 
cluding computation  techniques  and  material  modeling  methods.  Volume  II  out- 
lined a five-year  research  plan  intended  to  meet  the  objectives  of  the  non- 
destructive pavement  evaluation  procedure  outlined  above.  The  following  areas 
of  reasearch  were  proposed: 

(1)  Computer  programs  for  the  structural  analysis  of  pavements 

(2)  Engineering  studies  for  pavement  material  modeling 

(3)  Development  of  vibratory  test  equipment  and  data  interpretation 
techniques 

(4)  Data  collection  on  aircraft  characteristics  as  they  relate  to  the 
structural  response  of  pavements 

(5)  Test  and  evaluation  studies  to  validate  the  resulting  pavement  evalu- 
ation procedure 


RESEARCH  PROGRAM  (1970-1974) 
Computer  Code  Development 


The  principal  effort  on  the  computer  code  development  was  performed  at  the 
Civil  Engineering  Laboratory  (formerly  the  Naval  Civil  Engineering  Laboratory). 


1.  Nielsen,  J.  P.,  Rational  Pavement  Evaluation  - Review  of  Present  Technology , 

AFWL-TR-69-9,  Air  Force  Weapons  Laboratory,  Kirtland  Air  Force  Base,  New 
Mexico,  Vol . I,  October  1969;  Vol . II,  May  1970. 


The  resulting  computer  code  (AFPAV)  was  a nonlinear,  finite-element  program  cap- 
able of  treating  single  and  multiwheeled  landing  gears  (ref.  2).  The  code  is  ac- 
tually composed  of  three  main  programs--AFPRE  (data  input  and  element  generation), 

AFPAV  (equation  solving),  and  AFPOST  (recovery  of  results  for  plot  routines).  A 

detailed  AFPAV  user’s  manual  has  been  prepared  by  CERF  (ref.  3).  Recently,  a code 
known  as  AFCAN  (ref.  4)  was  prepared  to  greatly  reduce  the  manhours  and  potential 
errors  which  occur  when  AFPAV  is  implemented  in  its  original  form.  AFCAN  is  in- 
tended to  make  AFPAV  easy  to  use  in  routine  pavement  evaluation  studies.  The 
AFCAN  Code  also  contains  the  results  of  the  fourth  study  area  outlined  above. 

The  most  recent  modifications  of  the  AFPAV  Code  are  presented  in  section  3. 

Material  Modeling  Studies 

A particularly  potential  aspect  of  the  finite-element  technique  is  its  capabil- 
ity to  handle  generalized  constitutive  equations.  This  means,  for  example,  that 
the  nonlinear  stress/strain  behavior  of  subgrade  soils  can  be  considered  in  pave- 
ment evaluation  analyses.  Research  conducted  at  the  University  of  Kentucky  for 
AFWL  (ref.  5)  was  concerned  with  the  development  of  constitutive  equations  for 
base  course  and  subgrade  materials.  The  resulting  constitutive  equation,  known 
as  the  Hardin  hyberholic  shear/strain  equation , has  been  incorporated  into  the 
AFPAV  Code.  With  this  single  equation  it  is  possible  to  define  the  shear/strain 
relationship  for  base  course  or  subgrade  materials  for  a wide  range  of  condi- 
tions, including  soil  density,  saturation,  and  plasticity,  and  such  load  param- 
eters as  loading  rate  and  number  of  loading  cycles.  CERF  has  rewritten  the  gen- 
eral form  of  Hardin's  work  so  that  it  is  simply  inputted  to  the  AFPAV  Code.  This 
is  accomplished  by  adding  simple  soil  parameters  on  the  AFCAN  data  cards. 


2.  Crawford,  J.  E.,  An  Analytical  Model  for  Airfield  Pavement  Analysis , 

AFWL-TR-71 -70,  Air  Force  Weapons  Laboratory,  Kirtland  Air  Force  Base, 

New  Mexico,  May  1972. 

3.  Nielsen,  J.  P.,  AFPAV  Computer  Code  for  Structural  Analysis  of  Airfield 

Pavements , AFWL-TR-75-1 51 , Air  Force  Weapons  Laboratory,  Kirtland  Air 
Force  Base,  New  Mexico,  October  1975. 

4.  Crawford,  J.  E.,  Software  for  Everyday  Usage  of  AFPAV,  Technical  Memorandum 
M-51-76-06,  Civil  Engineering  Laboratory,  Port  Hueneme,  California,  March 
1976. 

5.  Hardin,  B.  0.,  Constitutive  Relationships  for  Airfield  Subgrade  and  Base 

Course  Materials,  Technical  Report  UKY  32-71-CE5,  College  of  Engineering, 
University  of  Kentucky,  Lexington,  Kentucky. 


Vibratory  Equipment  Development 


The  equipment  used  in  the  nondestructive  pavement  evaluation  procedure  (fig.  1) 
was  developed  at  CERF.  The  NDPT  van  is  8 ft  wide  and  35  ft  long  and  is  divided 
into  three  compartments  which  contain  all  the  test  equipment  necessary  for  non- 
destructive testing  of  pavements.  The  forward  generator  compartment  contains 
a 100-kW  diesel  generator  which  supplies  110/220-V  electrical  energy  for  the 
test  equipment.  The  center  compartment  contains  the  electromagnetic  vibrator, 
the  vibrator  oil-cooling  system,  the  power  supply  for  the  vibrator  field  coil, 
a step-up  transformer  to  convert  220  V to  440  V,  and  a hydraulic  system  for 
raising  and  lowering  the  vibrator. 

The  vibrator  is  lowered  by  two  hydraulic  jacks  into  position  on  the  pavement 
through  a hatch  in  the  floor  of  the  van.  A large  mass  is  used  to  keep  the  vi- 
brator baseplate  in  contact  with  the  pavement.  The  mass/vibrator  system  weighs 
6750  lb  and  comes  into  contact  with  the  pavement  through  a 12-in-diameter 
baseplate.  The  vibrator  consists  of  armature  coils  and  field  coils.  The  elec- 
trical field  moves  the  armature,  which  is  connected  thr  ugh  three  load  cells  to 
the  baseplate.  The  load  cells  are  used  to  measure  the  dynamic  load,  which  can 
be  varied  sinusoidally  about  the  static  load  up  to  5000  lb,  peak  to  peak.  A 
sweep  oscillator  servomechanism  is  used  to  continuously  vary  the  frequency  of 
the  dynamic  load  from  10  to  3500  Hz.  The  output  from  a velocity  pickup  located 
in  the  baseplate  is  integrated  to  obtain  the  displacement  of  the  baseplate. 


The  aft  compartment  of  the  van  contains  the  instrumentation,  recording  and 
monitoring  equipment,  and  the  power  amplifier  console  and  sweep  oscillator 
servo,  which  control  the  vibrator.  Instrumentation  is  provided  for  measuring 
the  force,  frequency,  acceleration,  velocity,  and  displacment  of  the  vibrator 
baseplate.  The  vertical  acceleration  of  the  pavement  at  selected  distances 
from  the  baseplate  can  be  measured  by  accelerometers  which  are  epoxied  to  the 
pavement  and  the  phase  angle  between  any  two  accelerometers  can  be  determined 
by  inputting  the  signals  to  a phase  computer.  Geophones  can  be  used  to  obtain 
data  on  the  deflection  basin.  The  recording  and  monitoring  equipment  consists 
of  a 14-channel,  FM  magnetic  tape  recorder/reproducer,  an  X-Y  recorder,  an 
oscilloscope,  a digital  voltmeter,  and  a digital  frequency  counter.  A servo- 
mechanism on  the  sweep  oscillator  is  used  to  hold  the  load  or  frequency  at  a 
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generator 

VIBRATOR  COMPARTMENT 

INSTRUMENTATION  COMPARTMENT  , 


Figure  1.  Nondestructive  Pavement  Test  (NDPT)  Van 


desired  level.  The  power  amplifier  console  contains  readout  units  for  load, 
displacement,  and  acceleration,  and  panel  meters  to  observe  any  malfunctions 
in  the  power  amplifier  unit.  The  compartment  is  also  equipped  with  a heating/ 
cooling  system  for  temperature  control.  Further  details  on  the  van  and  its 
operation  can  be  found  in  reference  6. 

Four  different  types  of  tests  can  be  performed  with  the  equipment  in  the  van. 
These  tests  are  as  follows: 

(1)  Load  Deflection  Test  at  Constant  Frequency 

Load  deflection  tests  can  be  performed  at  constant  frequency, 
usually  at  either  15  or  25  Hz.  The  load  is  gradually  increased 
to  the  limit  (5000  lb).  The  output  from  the  velocity  pickup 
located  in  the  baseplate  of  the  vibrator  is  integrated  to  obtain 
displacement.  Data  from  this  test  are  used  to  calculate  the 
dynamic  stiffness  modulus  of  e pavement. 

(2)  Resonant  Frequency  Test 

In  a resonant  frequency  test,  the  load  is  held  at  a constant, 
preselected  value  and  the  frequency  is  varied  to  obtain  a dis- 
placement/frequency plot.  This  test  provides  a response  curve 
which  can  be  used  to  determine  the  resonance  characteristics  of 
the  vibrator/ pavement  system. 

(3)  Deflection-Basin  Test 

In  a def lection-basin  test,  the  velocity  pickups  are  fixed  to 
the  pavement  on  a radial  line  through  the  center  of  the  load 
plate.  The  frequency  is  preselected  and  the  load  is  increased 
from  0 to  5000  lb.  The  velocity  outputs  are  recorded  on  mag- 
netic tape  and  later  played  back  through  the  integrator  to  ob- 
tain the  deflection  basin  for  any  applied  load  up  to  5000  lb. 


6.  Baird,  Glenn  T.,  et  al..  Instruction  Manual,  for  Mobile  Nondestructive 
Vibratory  Test  Equipment,  AFWL-TR-74-301 , Air  Force  Weapons  Laboratory, 
Kirtland  Air  Force  Base,  New  Mexico,  August  1975. 
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(4)  Wave  Propagation  Test 

The  wave  propagation  test  is  the  principal  test  used  in  the  Air 
Force  pavement  evaluation  procedure.  In  this  test  the  load  is 
held  constant,  usually  at  1000  lb,  and  a frequency  sweep  is  con- 
ducted from  10  to  3500  Hz.  The  signals  from  the  accelerometers 
are  recorded  and  later  processed  through  a phase  computer  to 
produce  a phase  angle/ frequency  plot , which  is  used  to  construct 
a dispersion  curve  for  each  test  station. 

See  section  4 for  details  on  data  reduction. 

Test  and  Evaluation  Studies 

Test  and  evaluation  studies  to  verify  the  NDPT  procedure  proceeded  slowly. 

The  principal  reason  for  this  cautious  approach  concerned  the  use  of  a concept 
(vibratory  testing)  which  had  the  theoretical  potential  of  being  useful  but 
had  only  been  used  in  a limited  manner  in  pavement  studies  and  often  with 
either  conflicting  or  noninterpretable  results.  An  initial  evaluation  effort 
(ref.  7)  was  concerned  with  the  construction  of  a special  test  pavement  at  the 
CERF  facility.  Three  test  sections  were  built--one  of  Portland  cement  con- 
crete and  two  of  asphalt  concrete  over  a deep  uniform  silty  clay  subgrade  with 
a CBR  in  the  range  of  8 to  12.  These  test  sections  served  as  a facility  on 
which  vibratory  testing  procedures  with  the  NDPT  van  were  formalized.  Phase 
velocity/wavelength  relationships  were  determined  from  the  vibratory  data  and 
these  were  used  to  construct  dispersion  curves.  These  tests  indicated  the 
usefulness  of  the  vibratory  procedure  and  suggested  certain  modifications  in 
equipment  and  data  interpretation  techniques. 

Subsequently,  the  test  van  was  taken  to  seven  airbases  in  the  southwest  for 
test  and  evaluation  under  field  conditions.  A series  of  CERF  letter  reports 
summarizes  these  studies.  This  road  tour  revealed  certain  equipment  problems-- 
mainly  the  inability  of  the  bench  stock  electronic  components  to  withstand 
road  vibrations.  On  the  positive  side  the  tour  verified  the  suggested  test 
procedures  and  revealed  that  interpretable  dispersion  curves  could  be  obtained. 

7.  Rao,  H.  A.  B.,  Nondestructive  Evaluation  of  Airfield  Pavements  (Phase  71, 

AFWL-TR-71 -75,  Air  Force  Weapons  Laboratory,  Kirtland  Air  Force  Base, 

New  Mexico,  December  1971. 
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Two  other  field  test  programs  which  provided  useful  data  later  in  the  develop- 
ment of  the  NDPT  program  were  the  pavement  distress  investigation  (ref.  8)  and 
the  WES  multiple-wheel,  heavy  gear-load  pavement  tests  (ref.  9).  The  major 
objective  of  the  former  was  to  establish  parameters  to  predict  remaining  life 
in  an  existing  pavement;  the  latter  provided  subgrade  failure  criteria  for  the 
wide-bodied,  jumbo  aircraft. 


PAVEMENT  RESEARCH  (1974  TO  PRESENT) 

Program  Evaluation  (July  to  October  1974) 

Before  July  1974,  the  various  elements  of  the  pavement  evaluation  development 
program  proceeded  primarily  as  independent  tasks.  This  was  necessary  because 
each  of  the  five  program  areas  called  for  research  and  development  which  pushed 
the  state-of-the-art  and  thus  called  for  creative  and  original  technological 
developments.  With  the  apparent  successful  development  of  the  various  program 
elements  by  July  1974,  efforts  were  directed  toward  the  organization  of  these 
elements  into  a composite  pavement  evaluation  procedure.  This  program  direc- 
tion occurred  hearly  simultaneously  with  a complete  change  in  CERF  technical 
personnel  in  the  soils  and  pavements  area. 

For  the  new  CERF  personnel  to  comply  with  the  directive,  it  was  necessary  for 
them  to  thoroughly  review  the  previous  efforts.  Initial  work  was  directed  to- 
ward a review  and  evaluation  of  some  of  the  NDPT  field  data.  This  was  a par- 
ticularly important  task  because  conflicting  results  had  been  informally  pre- 
sented and  there  was  some  concern  as  to  whether  or  not  vibratory  testing  could 
effectively  be  used  for  pavement  load  evaluation  studies. 


8.  O'Brien,  Ken,  and  Associates,  Distress  Criteria  for  Pavement.  Systems , 

AFWL-TR-73-226,  Air  Force  Weapons  Laboratory,  Kirtland  Air  Force  Base, 
New  Mexico,  April  1974. 

9.  Hammitt,  G.  M.,  et  al.,  Multiple-Wheel  Heavy  Gear  Load  Pavement  Testa, 

AFWL-TR-70-113,  Vol . IV,  Air  Force  Weapons  Laboratory,  Kirtland  Air 
Force  Base,  New  Mexico,  November  1971. 
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Another  area  of  concern  related  to  the  shear  moduli  obtained  from  vibra- 
tory testing  on  pavements;  in  general  the  modulus  values  were  much  too 
large.  After  a review  of  these  and  other  problem  areas,  CERF  presented 
a briefing  at  AFWL  late  in  August  1974  to  discuss  some  of  the  project 
findings.  The  following  significant  points  were  presented  during  this 
briefing: 

(1)  A review  of  test  data  collected  at  Altus  Air  Force  Base  indi- 
cated that  the  data  could  be  organized  to  conform  to  that  of 
a forced  vibrator  with  viscous  damping. 

(2)  Theoretical  considerations  revealed  that  NDPT  data  should  be 
corrected  before  shear  moduli  are  calculated.  (This  would 
yield  lower  moduli;  but,  the  form  of  the  correction  factors 
was  not  known  at  that  time.) 

(3)  Differences  in  the  NDPT  data  collected  with  the  CERF  van 
and  that  collected  by  WES  (ref.  10)  could  be  explained  on 
the  basis  of  certain  theoretical  considerations  presented 
by  Lysmer  (ref.  11). 

As  a result  of  the  above  points,  some  of  the  former  concerns  relative  to  NDPT 
testing  were  eliminated  and  CERF  took  the  position  that  sufficient  data  were 
available  to  form  a provisional  pavement  evaluation  procedure.  The  proce- 
dure recommended  was  a structural  approach  in  which  nondestructive  pavement 
testing  would  provide  in-situ  material  properties  data.  These  data  and 
the  aircraft  load  characteristics  would  be  inputted  to  the  AFPAV  Code  for 


10. 


11 


Hall,  J.  W.  and  Green,  J.  L.,  Nondestructive  Vibratory  Tea  tin; j of  Air- 
port Pavements,  Vol.  I,  Evaluation  Methodology  and  Experimental  Test 
Results , FAA  Report,  U.S.  Department  of  Transportation,  Federal 
Aviation  Administration,  Washington,  D.C.,  December  1973. 

Lysmer,  J.,  Vertical  '-lotion  of  Rigid  Footings , Contract  Report  No. 
3-115,  University  of  Michigan,  ORA  Project  05366,  Ann  Arbor,  Michiqan, 
June  1975. 
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structural  analysis,  and  then  certain  critical  stresses  and  strains  would 
be  reviewed  with  respect  to  material  failure  properties  to  predict  permitted 
aircraft  coverages.  Also  in  the  briefing,*  certain  problem  areas  that  needed 
additional  review  were  outlined  and  an  approach  to  meet  the  directive  of  or- 
ganizing the  program  elements  into  a workable  pavement  evaluation  procedure 
was  suggested.  Following  this  briefing  CERF  was  directed  to  formulate  a 
specific  plan  to  accomplish  this  objective. 

NDPT  Program  (October  1974  to  November  1975) 

In  October  1974  CERF  negotiated  a work  unit  with  AFWL  to  continue  the  devel- 
opment of  the  NDPT  project.  The  thrust  of  this  effort  was  directed  toward 
organizing  the  provisional  pavement  evaluation  procedure  suggested  as  a re- 
sult. of  the  previous  three  months  of  work  on  the  project.  This  new  work 
unit  directed  CERF  to  thoroughly  evaluate  the  electronic  components  in  the 
NDPT  van  and  to  replace  with  more  reliable  components  those  equipment  items 
which  had  often  failed  on  road  trips.  An  additional  effort  was  concerned 
with  the  detailed  evaluation  of  the  NDPT  data  collected  on  the  earlier  eval- 
uation trips  (circa  1973).  This  review  indicated  that  some  of  the  earlier 
data  were  poor  because  of  equipment  malfunctions.  More  importantly,  however, 
a review  of  the  good  data  suggested  a procedure  by  which  the  shear  moduli 
could  be  corrected.  CERF  also  demonstrated  that  the  NDPT  data  could  be  used 
in  conjunction  with  the  AFPAV  Code  to  predict  pavement  deflections  which 
were  in  substantial  agreement  with  those  measured  on  the  test  runways  (ref. 
8).  CERF  presented  these  findings  in  a briefing  report  to  AFWL  on  January  2, 
1975  (ref.  12).  Although  these  results  were  encouraging,  CERF  suggested 
that  a second  field  evaluation  trip  would  verify  the  data  reduction  scheme 
and  provide  an  opportunity  to  test  and  evaluate  a wide  range  of  pavement 


Ir  ‘ 

This  briefing  was  also  presented  to  the  Commander,  AFCEC  and  the  Di- 
rectorate of  Engineering  Materials,  Tyndall  Air  Force  Base  on  October 
17,  1974  because  of  the  planned  transfer  of  AFWL/DEZ  activities  to 
AFCEC. 

12.  Nielsen,  J.  P.,  firth'!  of  > Path 'turn!  /•>><  limit  ion  Si/olein,  Proqress  Report 
No.  3 (T.D.  5.11/00),  Civil  Engineering  Research  Facility,  University 
of  New  Mexico,  Albuquerque,  Mew  Mexico,  December  30,  1974. 


types.  For  example,  the  previous  surveys  did  not  include  pavements  with  sta- 
bilized base  courses,  complex  pavement  cross-sections  involving  multiple  layers, 
or  those  resulting  from  overlay  operations.  Accordingly,  CERF  submitted  a test 
plan  (ref.  13)  to  AFWL  to  perform  pavement  evaluation  studies  at  six  Air  Force 
Bases:  Williams,  Luke,  Holloman,  Laughlin,  Cannon,  and  Kirtland.  This  plan 

was  approved  on  February  23,  1975;  the  field  work  was  completed  by  the  end  of 
May,  1975.  The  analyses  of  these  data  were  completed  in  November,  at  which 
time  the  authors  presented  a paper  (ref.  14)  which  summarized  the  results  of 
the  field  tests  at  the  dympov i urn  on  Hander,  trui^tive.  Toot  and  Evaluation  of  A ir- 
port  Pavemcnto  held  at  WES.  An  enlarged  discussion  of  these  results  was  pub- 
lished in  1976  (ref.  15).  These  tests  demonstrated  an  excellent  corrplation 
between  field-measured  deflections  (under  actual  aircraft  loading)  and  those 
predicted  by  the  AFPAV  Code  with  NDPT  data.  These  results  suggested  that  the 
data-interpretation  and  correction  procedures  were  valid  and,  thus,  a reason- 
ably valid  pavement  evaluation  procedure  could  be  formulated.  As  a result  of 
the  successful  evaluation  of  the  NDPT  data  from  the  second  field  trip,  plans 
were  made  to  transfer  the  NDPT  van  to  AFCEC.  Accordingly,  AFCEC  personnel  were 
given  training  on  the  use  of  the  van  and  the  van  was  officially  transferred  to 
AFCEC  on  October  9,  1975.  AFCEC  is  currently  using  the  NDPT  van  in  its  pave- 
ment evaluation  program. 

One  area  of  particular  concern  was  the  reduction  of  the  NDPT  data.  During 
the  field  testing  phase  of  a pavement  evaluation  study,  a large  amount  of 
data  is  easily  acquired.  The  reduction  of  these  data  represents  a signifi- 
cant investment  of  manhours,  so  much  so  that  CERF  was  directed  to  develop 


13.  Nielsen,  J.  P.,  CERF  letter  to  Major  George  D.  Ballentine,  AFWL/DEZ, 
Kirtland  Air  Force  Base,  New  Mexico,  February  20,  1975. 

14.  Nielsen,  J.  P.,  and  Baird,  Glenn  T.,  Air  Faroe  iv/rtem  far  Nondestruai  i.ve 
Tmtin.ii  of  Pavementn , Symposium  on  Nondestructive  Test  and  Evaluation  of 
Airport  Pavements,  Waterways  Experiment  Station,  Corps  of  Engineers, 
Vicksburg,  Mississippi,  November  1975. 

15.  Nielsen,  J.  P.,  and  Baird,  Glenn  T. , Honda n true t.i.ve  lavement  Load  Ha  Id  raj , 

Thirteenth  Paving  Conference,  University  of  New  Mexico,  Albuquerque, 

New  Mexico,  January  1976. 


computer  routines  for  automatic  reduction  of  the  NDPT  data.  This  problem 
was  thoroughly  reviewed  with  the  AnaLog-to-Uigital  Section  of  AFWL.  It  soon 
became  evident  that  although  the  field  data  could  be  digitized,  no  automatic 
(computer-based)  techniques  existed  or  could  be  developed  to  provide  the  engi- 
neering judgment  which  is  necessary  to  interpret  the  field  data  prior  to  com- 
putation by  the  AFPAV  Code.  Therefore,  CERF  developed  a laboratory-based, 
data-reduction  electronic  package  which  can  be  used  to  reduce  the  field  data 
independent  of  the  NDPT  van.  This  setup  permits  the  field  crews  to  gather 
NDPT  data  with  the  van,  while  a laboratory-based  crew  proceeds  with  the  prep- 
aration of  phase  angle/frequency  plots  used  in  the  evaluation  studies.  Appen- 
dix A contains  operational  instructions  for  the  data  reduction  package. 

NDPT  Developments  Since  November  1975 

The  thrust  of  the  work  subsequent  to  the  developments  through  November  1975 
has  been  concerned  with  the  finalization  of  the  NDPT  test  and  data  evaluation 
procedures  and  the  orderly  transfer  of  equipment  and  technology  to  AFCEC. 


SECTION  3 

AFPAV  CODE  MODIFICATIONS 

INTRODUCTION 

The  AFPAV  Code,  which  performs  stress/strain  calculations  on  pavements  subjec- 
ted to  aircraft  loading  (ref.  16),  has  been  used  extensively  by  CERF  during  the 
development  and  evaluation  of  the  NDPT  procedure.  AFPAV  is  a nonlinear,  finite- 
element  program  in  which  the  pavement  is  represented  as  a two-dimensional  ex- 
tended solid.  The  code  is  capable  of  treating  both  single  and  multiwheeled 
landing  gears;  tire  loadings  are  represented  by  a Fourier  series.  A prepro- 
cessor code,  AFPRE,  formats  pavement  and  load  data  for  AFPAV.  Unfortunately, 
however,  preparation  of  the  AFPRE  data. is  a very  cumbersome,  error-prone  task. 
Therefore,  to  make  the  code  easier  to  use  and  to  provide  a means  of  predicting 
the  fatigue  life  of  pavements,  two  major  modifications  were  added  to  the  AFPAV 
Code--a  preprocessor  code  known  as  AFCAN  and  material  fatigue  subroutines.  The 
latter  were  placed  in  subroutine  RESULT. 

The  computer  program  AFCAN,  written  by  John  Crawford  of  CEL  provides  a simpler 
means  of  obtaining  the  input  data  required  for  AFPRE  and  AFPAV.  As  originally 
written,  AFCAN  was  used  to  generate  the  AFPRE  data  cards  which  were  subsequent- 
ly used  to  implement  AFPAV.  AFCAN  also  contained  a data  check  subroutine  which 
insured  that  reasonable  and  complete  data  were  introduced  into  AFPRE.  However, 
before  AFCAN  could  be  used  for  operational  studies,  the  values  assigned  to  its 
variables  had  to  be  reviewed  to  assure  the  user  that  realistic  stress  analyses 
were  being  performed. 

The  code  modifications  described  in  this  report  include  the  enlargement  of  AFCAN 
to  cover  more  aircraft  than  it  originally  covered  and  the  incorporation  of  AFCAN 
into  the  code  so  that  AFPAV  would  operate  efficiently  as  a single  code.  Thus, 
AFCAN  permits  the  highly  sophisticated  AFPAV  Code  to  be  used  for  routine  pave- 
ment evaluation  work  when  a large  number  of  AFPAV  runs  must  be  made  for  a wide 
range  of  military  aircraft. 

16.  Crawford,  J.,  and  Pichumani,  R.,  Finite-Element  Analiftiin  of  Pavement 

Struetures  ’teinq  AFPAV  Code.  (Nonlinear"  Elastie  Ana  levin ) , AFWL-TR-74-71  , 

Air  Force  Weapons  Laboratory,  Kirtland  Air  Force  Base,  New  Mexico,  April 
1975. 
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Any  pavement  evaluation  program  includes  both  pavement  load  rating  and  an  esti 
mate  of  the  fatigue  life  of  the  pavement.  AFPAV  performed  the  former,  but  the 
latter  had  to  be  added  to  the  code.  The  material  fatigue  considerations  used 
in  this  addition  to  AFPAV  are  also  described  in  this  section. 


AFCAN  CODE 
General 

The  AFCAN  Code  increases  efficiency  in  the  data  preparation  required  to  imple- 
ment AFPRE  and  AFPAV.  AFCAN  has  a series  of  subroutines  which  contain  the  con 
stants  for  a particular  aircraft:  tire  print  data,  aircraft  load  characteris- 

tics, Fourier  series  parameters,  and  finite-element  data.  Other  AFCAN  subrou- 
tines generate  all  of  the  data  cards  relative  to  node  and  element  locations 
and  print  statements.  Although  AFCAN  significantly  decreases  the  engineering 
effort  required  to  implement  AFPRE,  it  does  so  by  assigning  constant  values  to 
the  many  variables  in  AFCAN,  to  the  mesh  geometry  variables  subsequently  used 
in  AFPRE,  and  to  the  locations  (elements)  for  which  stress  data  are  printed. 

To  implement  AFCAN,  the  user  prepares  a series  of  cards  (one  for  each  layer 
of  the  pavement  system)  which  indicate  the  thickness  and  material  property 
data  and  an  identification  card  which  calls  a particular  aircraft  subroutine. 
These  few  cards  constitute  the  entire  input  to  AFCAN,  which  generates  data 
(cards  or  tape)  for  a complete  run  through  AFPAV.  Since  the  AFCAN  data  are 
easily  prepared  by  a technician,  decisions  relative  to  finite  elements  and 
nonlinear  material  property  considerations  are  eliminated. 

Finite-Element  Mesh 

Details  relative  to  the  AFCAN  and  AFPAV  Codes  are  fully  documented  in  ref- 
erences 3,  4,  and  17.  The  data  preset  in  the  AFCAN  subroutines  consist  of 


17.  Crawford,  J.  E.,  Unpublished  user's  manual  on  AFCAN.  Contract  re- 
port to  be  prepared  for  AFCEC  by  CEL.  Publication  date  unknown. 


mesh  generation  instructions  and  aircraft  loading  characteristics,  including 
data  for  the  Fourier  series.  Each  aircraft  listed  in  AFCAN  has  its  own  sub- 
routine, which  contains  the  following  code  variables  (constant  for  each  air- 
craft): 

NUMCOS  = number  of  Fourier  series  cosine  terms 
NELOAD  = number  of  elements  under  the  tire  (x-y  plane) 

NUMELX  = number  of  elements  along  the  x-axis 

WLOAD  = wheel  load 

ZLEN  = half-period  of  the  Fourier  series 

These  are  fixed  variables;  others  such  as  layer  thickness,  material  proper- 
ties, and  nonlinear  material  parameters  are  introduced  on  the  AFCAN  data 
cards. 

The  AFPAV  Code  treats  the  pavement  as  a solid  with  the  surface  of  the  pavement 
lying  in  the  x-z  plane;  the  x-y  plane  contains  the  cross-section  of  the  pave- 
ment. The  thickness  of  each  layer  is  read  directly  into  AFCAN  from  cards  pre- 
pared by  the  user.  The  depth  of  the  mesh  is  set  at  144  in.  This  depth  was 
arbitrarily  selected,  and  it  can  easily  be  changed  if  necessary.  The  half- 
width of  the  mesh  (XLEN)  was  also  arbitrarily  established.  For  a single  wheel 
load,  XLEN  is  set  at  25  times  the  half-width  of  the  tire  (XLEND).  The  number 
of  columns  of  elements  (NUMELX)  is  set  in  each  aircraft  subroutine  in  AFCAN. 

The  selection  of  values  for  these  variables,  which  are  concerned  with  the  ge- 
ometry of  the  mesh,  is  based  on  engineering  judgment,  consideration  of  the  ac- 
curacy of  the  stress  results,  and  computer  costs.  In  general,  the  size  of  the 
elements  should  be  such  that  the  flexibility  of  the  pavement  is  properly  mod- 
eled. This  is  achieved  by  making  the  elements  small;  but,  this  results  in  a 
large  number  of  elements  and  computer  costs  which  are  excessive.  As  the  size 
of  each  element  is  increased,  the  system  becomes  less  flexible  and  the  comput- 
er costs  decrease.  No  theoretical  guidance  is  available  for  making  these  de- 
cisions. Computer  runs  on  typical  pavements  must  be  made  and  the  results  com- 
pared before  a satisfactory  mesh  is  achieved.  The  selection  of  values  for  mesh 
variables  is  a subjective  matter.  Values  which  produce  satisfactory  stress  re- 
sults at  reasonable  computer  fees  have  been  assigned  to  these  variables. 
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Fourier  Series  Load  Function 


In  AFPAV  the  tire  loads  are  applied  normal  to  the  x-z  plane.  Gear  assemblies 
are  situated  to  produce  symmetrical  loading  conditions  with  respect  to  the 
x-y  and  y-z  planes.  A Fourier  series  is  used  to  model  tire  loads  for  single 
and  multiwheeled  landing  gears.  However,  the  tire  loads  are  replaced  with  an 
equivalent  set  of  forces  which  are  applied  to  the  nodes  directly  under  each 
tire.  Because  the  gear  is  located  to  produce  symmetry,  the  Fourier  series 
takes  the  form  of  an  even  function  in  which  the  coefficients  are  given  by 

% ■ -f  r f(z>  'Ki1)2112 

Jo 

where 

f(z)  = tire  pressure  (PSINOM)  defined  over  the  tire  length  (ZLENLD) 
p = half-period  of  the  function  (ZLEN) 
n = number  of  terms  (NUMCOS) 

The  Fourier  series  coefficients  are  generated  in  AFPRE  for  prescribed  values 
of  NUMCOS  and  ZLENLD.  The  Fourier  series  functional  representation  of  the 
tire  pressure  takes  the  following  form: 

n=l 

where  Fu^  is  the  force  at  node  i which  produces  a displacement  in  the 
y-direction.* 

The  force  vector,  (FI,  contains  a nodal  force  for  each  degree-of-freedom  for 
each  node  for  each  harmonic  of  the  Fourier  series.  In  the  above  expression 
for  Fu^ , the  z-term  in  cos(niT/p)z  represents  the  z-station  (x-y  plane)  at 
which  results  are  to  be  calculated;  in  the  expression  for  the  Fourier  coef- 
ficient, z represents  the  interval  of  integration.  Note  for  this  latter 
expression  that 

T 

This  does  not  strictly  follow  Fourier  series  concepts.  NUMCOS  = 1 corres- 
ponds to  1 / 2aQ , the  first  term  of  the  series,  not  the  first  cosine  term. 
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f(z)  = PSINOM  if  0 < z < ZLENLD/2 
f ( z)  = 0 if  ZLENLD/2  < z < ZLEN 


Selection  of  Constants 


The  technique  used  to  select  AFCAN  constants  is  best  outlined  by  detailing  the 
procedure  for  a particular  aircraft.  However,  the  analysis  is  fully  applicable 
to  all  single-wheeled  aircraft.  Each  multiwheeled  aircraft  must  be  considered 
as  a separate  case  in  which  selection  of  the  AFCAN  constants  is  more  subjective 
than  in  the  single-wheeled  case.  Values  for  tire  pressure  (PSINOM),  wheel  load 
(WLOAD),  tire  length  (ZLENLD),  and  half-width  (XLENLD)  are  based  on  landing  gear 
data  for  the  particular  aircraft  at  maximum  take-of*  gross  weight  (ref.  13). 
These  are  included  in  AFCAN  as  the  standard  constants,  however,  other  values  for 
WLOAD  and  PSINOM  can  be  entered  into  AFCAN.  If  WLOAD  and  PSINOM  are  not  enter- 
ed, the  code  defaults  to  the  values  for  the  aircraft  at  basic;  miss-'  n weight. 

For  the  F-4  aircraft,  the  following  standard  constants  are  listed  or  calculated 
in  AFCAN: 

WLOAD  = 27,000  lb 
PSINOM  = 265  psi 
XLENLD  = 4.43  in 
ZLENLD  = 11.50  in 

If  other  values  of  WLOAD  and/or  PSINOM  are  entered  into  AFCAN,  new  values 
of  XLENLD  and  ZLENLD  will  be  calculated  by  the  code. 

Selection  of  the  variables  associated  with  the  Fourier  series  requires  care- 
ful study.  The  Fourier  series  is  used  so  that  multiwheeled  landing  gears  can 
be  considered  in  AFPAV  and  stress  calculations  in  the  x-y  plane  at  any  z- 
station  (z-coordinate)  can  be  determined.  To  accomplish  this,  the  wheel  loads 
are  considered  to  be  periodic  functions,  and  the  tire-pressure  diagrams  in  the 
y-z  plane  are  modeled  by  a Fourier  series.  As  indicated  above,  the  tire  pres- 
sure is  considered  to  be  uniform.  However,  any  functional  distribution  could 
be  modeled  with  the  AFPAV  Code;  this  is  a potential  of  the  AFPAV  Code  which 
has  not  been  investigated. 


18.  Hay,  D.  R.,  Aircraft  Characteristics  for  Pavement  Design  and  Evaluation , 

AFWL-TR-69-54,  Air  Force  Weapons  Laboratory,  Kirtland  Air  Force  Base, 

New  Mexico,  October  1969. 
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The  primary  variables  which  play  a dominate  role  in  the  tire-pressure  diagram 
modeling  are  NUMCOS  and  ZLEN--the  number  of  cosine  terms  in  the  Fourier  series 
and  the  half-period  of  tne  load  function,  respectively.  There  are  no  general 
guidelines  available  for  assigning  numerical  values  to  these  variables.  A 
suggested  maximum  value  for  ZLEN  can  be  obtained  by  considering  the  aircraft 
landing  gear  layout.  From  a theoretical  point  of  view  the  degree  to  which  the 
Fourier  series  fits  or  describes  the  tire-pressure  diagram  can  be  evaluated  by 
several  fitness  techniques.  However,  an  excessively  large  number  of  cosine 
terms  is  required  to  achieve  a mathematically  satisfactory  fit  because  the  duty 
cycle  (ratio  of  ZLENLD  to  ZLEN)  is  small.  Additionally,  computer  costs  become 
excessive  as  the  number  of  cosine  terms  increases.  To  solve  this  apparent  di- 
lemma, values  for  NUMCOS  and  ZLEN  were  selected  on  the  basis  of  the  response 
of  the  pavement  to  the  Fourier  load  function,  rather  than  by  a mathematical 
measure  of  fit.  Before  explaining  the  approach  used  to  select  values  for  these 
variables,  a brief  discussion  concerning  the  interaction  of  NUMCOS  and  ZLEN  is 
necessary. 

Figure  2 is  a plot  of  the  Fourier  load  series  for  ZLEN  = 15  x (ZLENLD/2)  and 
various  values  of  NUMCOS.  For  the  F-4  aircraft,  this  corresponds  to  a numer- 
ical value  for  ZLEN  of  86.22  in,  or  a period  of  172.44  in.  In  this  plot  the 
load  is  normalized  to  the  tire  pressure  (PSINOM).  If  NUMCOS  =1,  the  load 
function  corresponds  to  a step  function.  The  addition  of  more  cosine  terms 
increases  the  normalized  load  and  narrows  the  length  of  the  first  half-cycle. 
Ideally,  this  distance  should  be  half  the  tire  length  and  the  pulse  should  be 
rectangular  as  indicated  in  the  figure.  Note  that  negative  load  values  occur; 
these  indicate  that  tension  (uplift)  is  applied  to  the  pavement.  However,  in 
each  case  the  total  load  applied  to  the  pavement  is  the  same  (twice  the  inte- 
gral of  PSINOM  over  the  interval  0 to  ZLENLD/2).  With  a large  number  of  terms 
the  series  becomes  concave  about  the  axis  of  symmetry  because  of  the  addition 
of  the  higher  frequency  harmonics.  The  fit  is  good,  but  the  computer  costs 
are  excessive.  To  improve  this  situation,  a smaller  value  may  be  used  for 
ZLEN.  A review  of  the  equation  given  for  the  Fourier  coefficients  shows  that 
the  argument  of  the  cosine  varies  inversely  with  p(ZLEN).  The  term  n-n/p  is  the 
angular  velocity;  therefore,  the  frequency  is  n/2p,  or  n/2  x ZLEN.  Thus,  as 
ZLEN  decreases,  the  values  of  the  coefficients  of  the  Fourier  series  increase 
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as  long  as  niT/ p ( ZL ENLD/ 2 ) < tt/2.  This  is  shown  in  table  1.  Note  that  for 
ZLEN  = 86.22  in,  all  coefficients  are  positive.  This  simply  means  in  this 
case  that  for  the  first  15  terms,  n7T/p(ZLENLD/2)  < it.  For  ZLEN  = 28.74  in, 

15  terms  extend  the  angle  through  2n  radians  to  an  angle  in  the  neighbor- 
hood of  3it  radians  (hence  the  oscillatory  nature  of  the  Fourier  coefficients). 
However,  the  first  four  terms  of  this  series  are  larger  than  those  in  the 
former  series.  This  suggests  that  fewer  terms  would  be  suitable  for  the 
smaller  value  of  ZLEN;  figure  3 shows  the  load  function  plotted  for  NUMCOS 
= 5 and  the  data  in  table  1.  In  the  case  of  the  smaller  ZLEN,  the  normal- 
ized load  reaches  a higher  value  and  provides  a better  fit  to  the  tire  pres- 
sure. For  this  case,  a 15-term  Fourier  series  with  ZLEN  = 86.22  in  produces 
the  same  load  function  as  5 terms  and  ZLEN  = 28.74  in.  However,  caution 
must  be  exercised  in  the  interpretation  of  this  result.  ZLEN  cannot  be  made 
small  simply  by  reducing  the  number  of  terms  in  the  series;  it  must  be  remem- 
bered that  the  function  is  periodic.  Thus,  the  second  pulse  of  the  function 

must  be  sufficiently  removed  along  the  z-axis  so  as  to  have  little  or  no  in- 

) 

fluence  on  the  stress  calculations  in  the  neighborhood  of  the  wheel  load. 

To  select  values  for  the  variables  which  influence  the  Fourier  series  load 
function,  the  response  of  a pavement  to  the  load  function  was  reviewed.  The 
approach  used  was  to  select  a pavement  section  and  to  observe  how  certain 
computed  results  varied  as  a function  of  the  Fourier  variables  and  the  flex- 
ibility of  the  pavement.  Although  the  load  function  is  not  affected  by  the 
flexibility  of  the  pavement,  the  stress  computations  are,  and  it  is  necessary 
to  consider  the  interaction  between  these  so  as  to  obtain  satisfactory  stress 
computations  and  reasonable  computer  costs.  An  8-in-thick  Portland  cement  con- 
crete pavement  supported  by  an  8-in-thick  stabilized  base  and  a 10-in- thick 
subbase  over  a 5000-psi  (elastic  modulus)  subgrade  was  studied.  This  rather 
strong  pavement  was  chosen  because  of  its  extreme  rigidity;  if  this  pavement 
could  be  reasonably  modeled — that  is,  its  flexibility  properly  captured  in  the 
finite-element  mesh--other  softer  pavements  could  also  be  satisfactorily  mod- 
eled. Pavement  response  was  evaluated  in  terms  of  the  tensile  stress  and 
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Table  1.  Fourier  Coefficients  as  a Function  of  ZLEN 


Coefficient 

Number 


Hi] 


z 

LEN 

86.22  in 

28.74  in 

0.067 

0.200 

0.132 

0.374 

0.130 

0.303 

0.125 

0.202 

0.118 

0.093 

0.110 

0.  101 

-0.063 

0.091 

-0.086 

0.079 

0.063 

0.055 

0.043 


strain  in  the  wearing  course  and  the  normal  stress  and  compressive  strain  in 
the  subgrade.  These  parameters  were  selected  because  they  are  needed  in  the 
application  of  this  program--the  load  rating  of  pavements  from  , -'^destructive 
test  (wave  propagation)  data. 


Function  Versus  ZLEN 


Pavement  flexibility  is  influenced  by  the  number  of  elements  along  the  x-axis 
(NUMELX)  and  the  number  of  rows  of  elements  within  each  material  layer.  As 
explained  earlier  this  latter  consideration  is  accomplished  by  AFCAN  without 
any  input  from  the  user;  however,  the  variable  NUMELX  is  in  each  aircraft 
subroutine.  Figure  4 indicates  how  the  concrete  tensile  stresses  and  strains 
under  the  centerline  of  the  load  (F-4  aircraft)  vary  with  NUMELX  and  NUMCOS 
for  the  case  of  ZLEN  = 36.22  in.  The  stress  variations  are  easily  explained 
by  reviewing  the  characteristics  of  the  Fourier  load  function  as  NUMCOS  in- 
creases. Low  stresses  and  strains  correspond  to  low  normalized  loads;  these 
occur  when  a few  Fourier  terms  are  used  and  at  higher  numbers  of  terms.  Ten 
cosine  terms  yield  the  highest  stresses  for  the  variables  plotted  for  ZLEN  = 
36.22  in.  The  effect  of  increasing  NUMELX  from  10  to  18  is  an  approximate 
4-percent  increase  in  the  computed  tensile  stress  for  all  values  of  NUMCOS 
shown.  However,  as  shown  in  figure  5,  the  computer  costs  increase  signifi- 
cantly with  increasing  values  for  both  NUMELX  and  NUMCOS.  It  was  an  objec- 
tive of  this  study  to  organize  the  code  so  that  the  cost  per  run  per  air- 
craft would  not  exceed  $10.00  on  the  CDC-7600  computer.  NUMELX  was  set  at 
14  because  the  boundary  stresses,  a and  a , were  very  small  and  because 

x y 

the  bending-stress  diagram  within  the  concrete  layer  was  nearly  symmetrical. 
However,  for  NUMCOS  = 10  and  NUMELX  = 14  the  cost  per  run  was  nearly  $15.00. 
In  order  to  reduce  the  run  time,  the  problem  was  studied  further.  As  dis- 
cussed earlier,  the  Fourier  load  series  is  a function  of  the  variable  ZLEN. 
Figure  6 is  a plot  of  the  argument  of  the  Fourier  series  2 x ZLEN/ZLENLD 
versus  the  concrete  tensile  stress  and  strain.  As  discussed  above,  NUMELX 
= 14  appeared  to  sufficiently  model  the  flexibility  of  the  pavement.  Ac- 
cording to  figure  5 fewer  than  seven  cosine  terms  would  be  required  to  keep 
the  cost  below  $10.00  per  run  for  NUMELX  = 14.  Therefore,  NUMELX  = 14  and 
NUMCOS  = 5 were  used  (fig.  6).  This  figure  indicates  that  the  tensile  stress 
reaches  a maximum  when  the  argument  is  7.5.  For  the  F-4  aircraft  this  cor- 
responds to  ZLEN  = 43.11  in.  Figure  4 shows  that  although  this  lower  value 
of  ZLEN  produces  a 5-percent  increase  in  tensile  stress  and  a smaller  in- 
crease in  the  tensile  strain,  the  computer  costs  are  reduced  to  $5.50  per 
run.  Subgrade  normal  stresses  and  compressive  strain  are  affected  by  ZLEN, 
NUMELX,  and  NUMCOS  in  a manner  similar  to  that  for  concrete  tensile  stresses 
and  strains. 


Tensile  Strain,  in/in  x 10" 


1 


Figure  5.  CDC-7600  Computer  Costs  as  a Function  of  NUMELX  and  NUMCOS 


Analyses  of  the  type  outlined  above  were  made  for  all  aircraft  (both  single 
and  multiwheeled  landing  gears)  in  the  AFCAN  library.  As  indicated,  the 
values  assigned  to  the  variables  were  selected  to  provide  data  for  the  Air 
Force  nondestructive  pavement  test  and  load  evaluation  program  and  to  min- 
imize computer  costs.  Although  these  variables  produce  satisfactory  results 
from  the  structural  analysis  (AFPAV),  the  decisions  which  led  to  their  se- 
lection were  subjective  and  these  values  may  not  be  the  same  as  those  selec- 
ted by  another  engineer.  This  is  especially  true  with  multiwheeled  aircraft, 
since  larger  numbers  of  elements  are  required  to  adequately  model  the  aircraft/ 
pavement  interaction.  In  these  cases  trade-offs  between  pavement  response, 
load  modeling,  and  boundary  stresses  were  made  in  order  to  limit  computer  costs 
to  an  average  of  $25.00*  per  run,  except  for  the  C-5A  which,  even  for  conser- 
vative modeling,  costs  about  $50.00  per  run.  However,  at  this  time  further 
refinements  do  not  seem  appropriate,  at  least  until  better  pavement  fatigue 
criteria  become  available.  The  AFCAN  constants  assigned  to  each  aircraft  are 
tabulated  in  table  2. 

★ " 

The  cost  for  the  CDC-7600  computer  at  $610.00  per  system-hour. 
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Aircraft  WLOAD,  PSINOM,  NUMCOS  ZLEN , Cycles  per  Coverage 
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LOAD  RATING  SUBROUTINE 


The  AFPAV  Code  provides  for  a complete  structural  analysis  of  a pavement. 
However,  before  an  airfield  pavement  can  be  load  rated,*  known  material  fa- 
tigue characteristics  must  be  considered  along  with  the  results  of  the  struc- 
tural analysis.  The  capability  to  load  rate  a pavement  is  the  principal  ob- 
jective of  the  Air  Force  research  and  development  program  for  pavement  eval- 
uation. With  the  successful  development  of  the  NDPT  van  and  the  associated 
data-reduction  schemes,  a provisional  load-rating  technique  for  pavements 
was  formulated  and  placed  in  AFPAV. 


Failure  Criteria 


In  any  load-rating  procedure,  it  is  necessary  to  establish  the  most  probable 
failure  modes  for  the  pavement.  The  following  failure  criteria  were  selected: 

(1)  The  maximum  tensile  strain  in  the  bituminous  wearing  course 
should  not  exceed  a prescribed  level  of  strain  expressed  in 
micro  inches. 

(2)  The  maximum  tensile  stress  in  a Portland  cement  wearing  course 
should  be  below  the  modulus  of  rupture  of  the  concrete. 

(3)  The  vertical  compressive  strain  in  the  subgrade  should  not  ex- 
ceed a prescribed  level. 

To  implement  these  structural  procedures,  failure  criteria  for  bituminous 
surface  courses,  Portland  cement  concrete,  and  subgrade  materials  were  ex- 
tracted from  the  literature. 


Bituminous  Concrete.  --  The  fatigue  criteria  shown  in  figure  7 (ref.  19)  were 
used  in  the  development  of  a load-rating  procedure  for  bituminous  pavements. 
Although  several  other  criteria  are  available,  those  in  figure  7 are  useful 
in  that  the  tensile  strain  in  the  bituminous  layer  for  a particular  pavement 


* 

Load  rating  is  the  process  of  indicating  what  aircraft  and  how  many  annual 
movements  of  each  aircraft  can  safely  be  made  from  a given  airfield  pave- 
ment. 

19.  Monismith,  C.  L.,  and  McLean,  0.  B.,  "Structural  Design  Considerations," 
Pvooeediruiu , Association  of  Asphalt  Paving  Technologists,  Cleveland, 
Ohio,  1966. 
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and  aircraft  loading  can  be  determined  from  a structural  analysis.  The 
elastic  modulus,  E,  is  determined  from  the  nondestructive  field  test.  With 
these  data  it  then  becomes  a simple  matter  to  determine  from  figure  7 the 
number  of  load  repetitions  required  to  cause  failure  in  the  bituminous  wear- 
ing course  for  a particular  aircraft. 

Concrete.  --  The  Portland  Cement  Association  failure  criteria  for  Portland 
cement  concrete  (ref.  20)  were  adopted  for  concrete  pavements.  A graphical 
representation  of  the  concrete  failure  criteria  is  shown  in  figure  8. 

Subgrade.  --  Figure  9 is  a graphical  representation  of  the  fatigue  criteria 
selected  for  the  subgrade  materials.  These  relationships  were  taken  from 
data  presented  by  WES  from  its  test  section  studies  of  conventional  and 
heavy,  multiwheeled  aircraft,  from  Shell  Oil  Company  suggested  criteria, 
and  from  data  presented  by  The  Asphalt  Institute.  Three  fatigue  relation- 
ships have  been  used  to  cover  the  conditions  of  Portland  cement  concrete 
pavements  and  strong  bituminous  pavements  (E  > 350,000  psi),  weak  bitu- 
minous pavements  (E  < 200,000  psi),  and  heavy,  multiwheeled  aircraft 
(B52,  C5A,  etc.). 

Code  Revisions  for  Fatigue  Analysis 


To  incorporate  a pavement  fatigue  analysis  into  the  AFPAV  Code,  certain 
changes  and  additions  had  to  be  made  to  AFPAV.  To  accomplish  these  changes 
a new  code,  known  as  PREDICT,  was  prepared.  PREDICT  performs  structural 
analyses  on  many  problems  (up  to  20)  and  it  also  contains  a subroutine  to 
perform  mixed  traffic  analysis. 

PREDICT  is  a single  program  which  in  essence  is  composed  of  AFCAN,  AFPRE,  AFPAV, 
and  the  fatigue  and  mixed  traffic  subroutines.  Some  of  the  capabilities  pre- 
sent in  AFPAV  (ref.  3),  such  as  formating  of  data  for  RESP0N  (the  plot  pack- 
ages DEFLEC  and  CONTOUR),  have  not  been  included  in  PREDICT.  However, 


20.  Packard,  R.G.,  "Design  of  Concrete  Airport  Pavement,"  Engineering 
Bullet-in , Portland  Cement  Association,  1973. 
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Asphalt  Concrete  Tensile  Strain,  in/in 


Modulus  of  Elasticity,  psi 
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strain  below  70  pin/in 
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Number  of  Repetitions  to  Failure 


Figure  7.  Fatigue  Criteria  for  Bituminous  Pavements 
[after  Moni smith  (ref.  19)] 
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Figure  8.  Fatigue  Criteria  for  Concrete  Pavements 
[after  Packard  (ref.  20)] 


PREDICT  makes  maximum  use  of  previous  developments  associated  with  AFPAV 
and  AFCAN. 

All  elements  of  the  fatigue  analysis  have  been  placed  in  the  AFPAV  subroutine 
RESULT.  As  a first  step  in  this  analysis,  a subroutine  was  written  to  search 
for  the  maximum  stress  or  strain  values  induced  in  the  pavement  by  the  air- 
craft loading.  This  subroutine  searches  a series  of  prescribed  x-y  planes 
(z-stations)  for  the  critical  fatigue  factors.  The  results  of  this  search 
are  printed  as  part  of  the  computer  output  in  a format  which  indicates  the 
z-station  (ZPLANE)  at  which  the  values  are  maximum.  If  the  wearing  course 
is  a bituminous  surface,  the  maximum  tensile  strain  in  the  surface  course 
will  be  printed  as  well  as  the  maximum  compressive  strain  in  the  subgrade; 
if  the  pavement  has  a concrete  surface,  the  maximum  tensile  stress  in  the 
concrete  as  well  as  the  maximum  subgrade  compressive  strain  will  be  printed. 
These  stresses  and  strains  are  used  in  the  fatigue  analyses. 

A^orithms  were  written  for  the  fatigue  relationships  presented  in  figures 
7,  8,  and  9.  These  algorithms  were  arranged  into  three  subroutines  so  that 
the  number  of  cycles  to  failure  for  each  material  (surface  and  subgrade)  can 
be  calculated.  With  the  concept  of  cycles  per  coverage  (ref.  18),  the  num- 
ber of  operations  (OPS)  permitted  for  each  material  is  printed  following  the 
stress  and  strain  data.  An  operation  is  defined  as  a single  movement  of  the 
aircraft  (departure)  at  basic  mission  weight.  Obviously,  the  lower  of  the 
two  reported  OPS  controls  for  the  particular  aircraft  and  runway  under  con- 
sideration. The  controlling  pavement  component  (subgrade  or  wearing  course) 
is  also  indicated  on  the  printout.  This  information  is  useful  in  planning  a 
program  to  upgrade  the  load  capacity  of  the  pavement. 

Mixed  Traffic  Considerations 

Because  more  than  one  type  of  aircraft  will  most  certainly  use  a given  run- 
way or  taxiway,  the  influence  of  mixed  aircraft  usage  was  also  considered. 

The  procedure  developed  consists  of  establishing  a mixed  traffic  factor 
based  on  the  limiting  value  of  OPS,  the  aircraft  weight  and  the  cycles  per 
coverage  for  each  aircraft.  The  fixed  factors  for  this  analysis  have  been 


placed  in  the  aircraft  subroutines  in  AFCAN.  The  controlling  mixed  traffic 
factor  is  printed  in  the  computer  output  along  with  the  stress/strain  data 
and  the  OPS  information.  A typical  computer  summary  is  shown  below: 


SHAW  AFB  RW  22R  SOOTH  END 
F4 


OPS-CONCRETE  OPS-ASPHALT  OPS-SUBGRADE 

85107.  0.  206943. 

MIXED  TRAFFIC  FACTOR  FOR  CONCRETE  = .0000008782 

MIXED  TRAFFIC  FACTOR  FOR  SUBGRADE  = .0000003612 

CONTROLLING  OPERATIONS  FOR  F4  = CONCRETE 

CONTROLLING  MIXED  TRAFFIC  FACTOR  FOR  F4  = CONCRETE 


SHAW  AFB  RW  22R  CENTER 
F4 


OPS-CONCRETE  OPS-ASPHALT  OPS-SUBGRADE 

1.  0.  45692. 

MIXED  TRAFFIC  FACTOR  FOR  CONCRETE  = .0747384155 

MIXED  TRAFFIC  FACTOR  FOR  SUBGRADE  = .0000016357 

F4  CANNOT  USE  THIS  FACILITY 


SHAW  AFB  RW  22R  NORTH  END 
F4 


OPS-CONCRETE  OPS-ASPHALT  OPS-SUBGRADE 

1338000.  0.  6690000. 

MIXED  TRAFFIC  FACTOR  FOR  CONCRETE  = .0000000559 

MIXED  TRAFFIC  FACTOR  FOR  SUBGRADE  = .0000000112 

CONTROLLING  OPERATIONS  FOR  F4  = CONCRETE 

CONTROLLING  MIXED  TRAFFIC  FACTOR  FOR  F4  = CONCRETE 


To  evaluate  the  effect  of  mixed  traffic,  the  user  simply  multiplies  the 
mixed  traffic  factor  for  each  aircraft  by  the  expected  number  of  annual  op- 
erations. The  sum  of  all  these  multiplications  must  not  exceed  one.  Each 
aircraft  which  is  expected  to  use  the  runway  must  be  processed  through 
PREDICT. 


If  the  critical  stress  or  strain  induced  in  the  pavement  by  a particular  air- 
craft falls  below  the  endurance  limit  of  the  material,  the  number  of  repeti- 
tions to  failure  is  as  shown  in  table  3.  These  values  are  then  converted  to 
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operations.  If  the  pavement  is  overstressed,  a statement  which  indicates 
that  the  particular  aircraft  cannot  use  the  facility  is  printed;  if  the 
computed  number  of  operations  is  between  1 and  10,  a statement  which  indi 
cates  that  only  emergency  use  is  permitted  is  printed. 


Table  3.  Limiting  Failure  Criteria  for  Paving  Materials 


Factor 

Material 

Asphalt 

Concrete 

Subgrade 

Endurance  Limit 

70  x 10-6  in/ in 

54  % Modulus  of 

7 x lO”"  in/in 

Rupture 

Repetitions 

106 

10s 

5 x 105 

SECTION  4 

DATA  INTERPRETATION  PROCEDURE 


to  develop  a procedure  for  interpreting  phase  angle/ 
frequency  plots  and  resulting  dispersion  curves.  This  procedure  is  to  be 
used  by  AFCEC  personnel,  along  with  the  NDPT  van  and  PREDICT,  to  perform 
pavement  evaluation  studies. 

A wave-propagation  technique  to  determine  the  load-carrying  capability  of 
airfield  pavements  has  been  developed  at  CERF.  This  technique  consists  of 
placing  a vibrator  rn  the  surface  of  a pavement  and  sweeping  the  load  through 
a frequency  range  of  10  to  3500  Hz.  The  propagation  velocity,  which  varies 
with  frequency,  is  determined  from  measurements  of  the  phase  angle  taken 
during  the  frequency  sweep  between  accelerometers  fixed  to  the  surface  of 
the  pavement.  This  variation  of  stress  wave  velocity  with  frequency  is  termed 
dispersion . A dispersion  curve  shows  how  the  propagation  velocity  changes 
with  frequency.  The  curve  usually  is  made  up  of  various  line  segments,  the 
terminal  or  peak  points  of  which  represent  the  critical  velocities  of  the 
propagating  stress  waves.  These  velocities  are  used  to  determine  the  elastic 
properties  of  the  different  layers  of  the  pavement  system.  These  properties 
are  then  used  in  a multilayered,  nonli'near,  finite-element  computer  code  to 
determine  the  structural  response  of  the  pavement  to  various  aircraft  load- 
ings. The  results  of  this  analysis  along  with  pavement  fatigue  data  provide 
a means  of  predicting  the  remaining  service  life  of  the  pavement. 


INTRODUCTION 


A study  was  undertake 


PHASE  ANGLE/FREQUENCY  PLOTS 

Data  from  several  airbases  tested  with  the  NDPT  van  in  April  1975  were 
used  (ref.  21)  in  a study  of  phase  angle/ frequency  plots.  The  phase 

21.  Nielsen,  J.  P.,  Airfield  Tost  Plan,  CERF  Letter  Report,  Aerospace  Fac- 
ilities Branch,  Air  Force  Weapons  Laboratory,  Kirtland  Air  Force  Base, 
New  Mexico,  February  20,  1975. 
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angle/frequency  plots  were  constructed  by  the  laboratory-based,  data -reduct ion 
system  developed  at  CERF  for  AFCEC  (ref.  22).  This  reduction  system  is  essen- 
tially the  same  as  that  in  the  NDPT  van  (ref.  6),  except  that  a larger  X-Y 
recorder  is  used  to  plot  the  data. 

Phase  angle/frequency  plots  were  obtained  for  different  pairs  of  accelero- 
meters to  determine  variations  in  the  data  at  different  distances  from  the 
vibrator.  Accelerometers  1 and  2 were  located  side-by-side  at  a distance  of 

2 ft  from  the  center  of  the  vibrator;  3 and  4 were  similarly  placed  4 ft 

from  the  center  of  the  vibrator;  5 and  6 were  located  6 ft  from  the  center 
of  the  vibrator  (ref.  6).  The  accelerometers  are  positioned  in  pairs  for 
two  reasons:  (1)  for  backup;  i.e.,  in  case  an  accelerometer  should  fail, 

data  could  still  be  obtained  by  the  other,  and  (2)  for  calibration;  i.e., 
assuming  that  a symmetrical  or  spherical  plane  wave  is  created  by  the  vibra- 
tor, a phase  angle  of  0 deg  should  result  between  adjacent  accelerometers 
and  thus  indicate  that  both  accelerometers  are  responding  to  the  same  sur- 
face excitation. 

Calibration  curves  for  the  three  pairs  of  accelerometers  are  shown  in  figure 
10  for  the  low-frequency  sweep  (10  to  350  Hz)  and  the  high-frequency  sweep 
(10  to  3500  Hz)  on  the  Portland  cement  concrete  pavement  at  Cannon  Air  Force 
Base.  The  calibration  curves  fcr  the  low-frequency  sweep  show  a near  0-deg 
phase  shift  between  the  paired  accelerometers.  The  data  from  accelerometers 

3 and  4 show  a greater  deviation  from  the  0-deg  phase  shift  line  than  the 
data  from  accelerometers  1 and  2 or  5 and  6.  The  high-frequency  sweep  data 
show  that  the  best  calibration  is  between  accelerometers  1 and  2 and  the 
poorest  is  between  5 and  6.  This  trend  was  observed  on  most  of  the  data 
and  can  be  used  as  an  indicator  of  the  quality  of  the  data.  Most  isolated 
spikes  are  believed  to  be  pavement  and/or  equipment  resonant  points;  some 
simply  cannot  be  explained. 

Phase  angle/ frequency  plots  were  prepared  for  12  test  locations  along  runway 
03-21  at  Cannon  Air  Force  Base.  Accelerometer  pairs  1 and  3,  and  2 and  3; 


22.  Nielsen,  J.  P.,  Development  of  Pavement  Evaluation  System,  Progress  Re 
port  No.  2 (T.D.  5.01/00),  Civil  Engineering  Research  Facility,  Kirt- 
land  Air  Force  Base,  New  Mexico,  January  31,  1976. 


Low-Frequency  Sweep 

bration  Curves  on  Portland  Cement  Concrete  Pavement  (1  of  2) 
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Calibration  Curves  on  Portland  Cement  Concrete  Pavement 


3 and  5,  and  4 and  5;  and  1 and  5,  and  2 and  5 were  used.  Typical  low- 
frequency  plots  are  shown  in  figure  11;  figure  12  shows  the  high-frequency 
data.  After  considerable  review  it  was  decided  to  interpret  these  plots 
by  following  the  general  trend  of  the  phase  angle/frequency  data  by  fair- 
ing through  the  spikes  (noise).  The  final  interpretation  of  the  data  is 
shown  as  a dashed  line.  These  interpretations  are  based  on  an  averaging 
of  the  data;  the  abrupt  slopes  where  the  plotter  shifts  through  360  deg 
were  ignored.  Data  points  along  the  interpreted  curve  were  used  to  es- 
tablish the  dispersion  curves  shown  in  figures  13a,  b,  and  c.  (See 
section  5 for  full  details.) 


As  part  of  the  phase  angle/frequency  plot  interpretation,  it  was  observed  that 
the  high-frequency  (short  wavelength)  portion  of  the  dispersion  curves  was 
more  consistent  for  the  data  from  accelerometers  1 and  3 and  2 and  3.  More 
consistent  low-frequency  data,  which  represent  subgrade  conditions,  resulted 
from  accelerometers  1 and  5 und  2 and  5.  If  the  low-frequency  (25  to  125  Hz) 
data  froi  accelerometers  1 and  5 are  combined  with  the  high-frequency  (250  to 
3500  Hz)  data  from  accelerometers  1 and  3,  a satisfactory  dispersion  curve  can 
be  developed  (fig.  14).  Because  the  high-frequency  (short  wavelength)  vibra- 
tions are  essentially  contained  in  the  surface  and  near-surface  layers,  it 
appears  that  the  closer-spaced  accelerometers  (i.e.,  1 and  3 or  2 and  3)  pro- 
vide a more  accurate  phase  angle/frequency  relationship  for  these  material 
layers.  A parallel  rational  was  applied  to  the  low-frequency  data,  which  re- 
present the  deeper  material  layers.  Thus,  greater  spacing  is  required  be- 
tween accelerometers  (i.e.,  4 ft  for  accelerometers  1 and  5 or  2 and  5)  to 
obtain  good  subgrade  data.  This  is  similar  to  a conventional  seismic  field 
exploration  in  which  data  obtained  by  closely  spaced  geophones  yield  infor- 
mation concerning  the  shallow  materials  and  those  with  a greater  spacing  give 
data  on  the  deeper  materials  (ref.  23). 


23.  Dobrin,  Milton  B.,  Introduction  to  Geophysical  Prospecting,  2nd  Edition, 
McGraw-Hill,  New  York,  1960. 
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Figure  11.  Typical  Low-Frequency  Phase  Angle/Frequency  Plots  (2  of  3) 
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Frequency  Phase  Angle/Frequency  Plots  (3  of  3) 
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Figure  12.  Typical  High-Frequency  Phase  Angle/Frequency  Plots  (2  of  3) 
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Figure  13.  Typical  Dispersion  Curves  (2  of  3) 


0009 


Accelerometers 

Accelerometers 

Accelerometers 


o 

o 

o 

o 

o 

o 

o 

o 

o 

CO 

C\J 

' 

Wavelength,  ft 

Figure  14.  Composite  Dispersion  Curve 
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DISPERSION  CURVES 

Dispersion  curves  were  developed  in  the  above  manner  for  the  test  locations 
at  Cannon  Air  Force  Base  and  nine  test  locations  at  Laughlin  Air  Force  Base. 
(Laughlin  has  an  asphaltic  concrete  surface  with  stabilized  base  courses.) 

The  interpretative  procedure  proved  to  be  very  reliable  when  applied  to  Port- 
land cement  concrete  pavements;  more  data  scatter  occurred  for  the  asphalt 
concrete  pavements.  Subsequent  application  of  this  technique  to  data  collec- 
ted at  Shaw  Air  Force  Base  proved  the  technique  to  be  equally  reliable.  Thus, 
the  procedure  is  recommended  for  the  pavement  evaluation  program.  A complete 
NDPT  data  reduction  guide  is  presented  in  section  5. 

Although  interpretation  of  the  phase  angle/frequency  plots  requires  some  ex- 
perience, analysis  of  the  dispersion  curves  is  the  critical  part  of  the  entire 
data-reduction  system.  Therefore,  a study  was  initiated  to  determine  the  sen- 
sitivity of  the  dispersion  curves  to  changes  in  material  properties  and  layer 
thicknesses.  To  accomplish  this,  a computer  code,  known  as  the  BRITISH  Code, 
was  used  to  obtain  solutions  to  the  differential  equations  of  motion  which 
are  applicable  to  this  problem.  The  known  compressiona 1 and  shear  velocities, 
mass  density,  and  thickness  of  each  material  layer  are  inputted  to  the  code 
in  a normalized  form.  All  velocities,  densities,  and  layer  thicknesses  are 
normalized  to  those  of  the  surface  layer.  The  user  must  also  specify  the  fre- 
quency range,  a frequency  increment,  and  an  estimate  of  the  complex  root  to 
the  differential  equations.  The  estimated  root  of  the  equation  is  used  as  an 
initial  value  in  the  solution  procedure.  The  code  then  uses  an  iterative  pro- 
cedure to  solve  for  the  phase  velocity  at  each  frequency. 

Several  dispersion  curves  reported  by  Jones  (ref.  24)  were  duplicated  with 
the  BRITISH  Code  as  shown  in  figure  15.  Figure  15a  shows  the  dispersion 
curve  for  a two-layered  system  (i.e.,  a surface  layer  over  a semi -inf ini te 
medium),  and  figure  15b  shows  the  dispersion  curve  for  a three-layered  sys- 
tem (i.e.,  a medium-speed  (phase  velocity)  layer  over  a high-speed  layer 


24.  Jones,  R.,  Thrower,  E.  N.,  and  Gatfield,  E.  N.  "Surface  Wave  Method," 

Peoaeodinpo  of  the  2nd  Into1.  Con]',  on  the  HLwurtupii  / dei'ljn  of  Asphalt 

1‘ai'emr.ritf , University  of  Michigan,  August  7-11  , 1967. 
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Velocity/Shear  Velocity  (Layer  1) 


Flexural  Mode 
(Composite  of 
Layers  1 and  2) 


Wavelength/Thickness  of  Layer  1 


(b)  Three-Layered  System 


Layer 

_ V\ 

1 

2 

3 

1.8 

4.475 

0.6833 

Figure  15.  Duplicate  of  Dispersion  Curves  by  Jones  (2  of  2) 


over  a low-speed,  semi-infinite  medium).  These  dispersion  curves  were  obtained 
only  after  a considerable  trial-and-error  effort  to  satisfactorily  duplicate 
Jones'  curves.  The  dispersion  curves  were  highly  dependent  on  the  input  data, 
specifically  the  initial  estimates  of  the  roots  of  the  equation.  This  is  il- 
lustrated by  the  dispersion  curve  results  (figs.  15b  and  16).  Figure  16  is 
the  computed  dispersion  curve  using  data  identical  to  that  for  figure  15b  but 
with  a different  estimate  of  the  imaginary  root  for  the  initial  frequency. 

No  analytical  method  could  be  developed  by  which  the  estimate  of  the  input 
root  could  be  improved  so  as  to  eliminate  the  problem  revealed  by  these  two 
figures.  It  became  evident  that  consistent  dispersion  curves  could  not  be 
obtained  with  the  BRITISH  Code  in  its  present  form.  To  further  illustrate 
the  inconsistent  results  from  the  code,  the  dispersion  curve  for  a third  prob- 
lem involving  a high-speed  layer  over  a medium-speed  layer  over  a low-speed 
layer  is  presented  in  figure  17  (ref.  24).  Figure  18  shows  the  best  solution 
that  could  be  obtained  for  this  problem  with  the  BRITISH  Code.  The  surface 
layer  was  then  divided  into  four  layers  of  equal  thickness  with  identical  ma- 
terial properties;  the  second  layer  was  divided  into  two  layers  with  identical 
properties.  The  total  thickness  of  the  system  was  not  changed;  but,  the  three- 
layered system  was  modified  to  appear  as  a seven-layered  system.  The  input 
frequency  and  estimated  roots  were  the  same  as  those  used  to  construct  figure 
13  (the  three-layered  solution).  Figure  19  shows  the  dispersion  curve  for  the 
seven-layered  system.  A considerable  difference  can  be  seen  between  these  two 
plots. 

Based  on  the  foregoing  results,  it  was  decided  that  the  BRITISH  Code  could 
not  be  used  with  confidence  as  a guide  to  interpret  field  dispersion  curves. 
Other  researchers  (ref.  25)  who  have  had  considerable  experience  with  other 
computer  codes  for  calculating  dispersion  data  have  had  little  success  in 
applying  such  codes  to  a system  in  which  the  stiffness  decreases  with  depth 
(i.e.,  pavement).  It  appears  that  present  codes  are  theoretically  valid, 
but  discrepancies  result  when  these  codes  are  applied  to  a pavement  system. 

In  addition,  there  is  some  concern  about  the  theoretical  basis  of  these 


25.  Personal  communication  with  Professor  John  Lysmer,  Civil  Engineering  De- 
partment, University  of  California,  Berkeley,  California,  February  4,  1976. 
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Figure  16.  Alternate  Solution  for  Three-Layered  System  for 
Different  Initial  Estimate  of  Imaginary  Root 
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Figure  19.  Dispersion  Curves  for  Seven-Layered  Representation 
of  System  Considered  in  Figure  18 
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codes  as  it  relates  to  the  NDPT  method  developed  by  CER>".  A comparison  of  the 
successful  dispersion  curves  computed  by  the  BRITISH  Code  (figs.  15,  16,  18, 
and  19)  and  those  computed  from  the  NDPT  data  (figs.  13  and  14)  indicates  very 
different  responses  of  the  pavement.  The  BRITISH  Code  finds  solutions  to  the 
differential  equations  of  motion  for  the  different  modes  of  vibration  (i.e., 
longitudinal,  flexual)  for  a layered  system;  whereas,  the  NDPT  data  are  an 
evaluation  of  the  componite  response  of  the  pavement.  It  can  be  seen  from  the 
NDPT  dispersion  curve  (fig.  14)  that  only  one  curve  with  various  segments  ex- 
ists. This  indicates  that  the  pavement  system  possesses  only  one  real  mode  of 
vibration  (principally  a flexural  mode)  as  opposed  to  the  theoretical  consid- 
erations of  the  BRITISH  Code.  Considerable  work  is  required  to  improve  and 
develop  the  BRITISH  Code  and  other  existing  codes  so  as  to  accurately  model 
a pavement  in  order  to  study  dispersive  media  or  systems  with  decreasing 
rigidity. 


SECTION  5 

DATA-REQUCTION  TECHNIQUES 


PHASE  ANGLE/FREQUENCY  PLOTS 

1 Phase  angle/frequency  plots  must  be  obtained  with  either  the  NDPT  electronics 

in  the  van  or  the  laboratory,  data-reduction  instrumentation.  Plots  for  re- 
duction should  include  the  following: 

Plot  No.  Frequency  Range,  Hz  Accelerometer  Pairs 

1 10  to  350  1 and  3,  2 and  3,  1 and  2 

2 10  to  350  1 and  5,  2 and  5,  5 and  6 

3 10  to  3500  1 and  3,  2 and  3,  1 and  2 

4 10  to  3500  1 and  5,  2 and  5,  5 and  6 

In  addition,  plots  of  the  phase  angle  in  the  +180-deg  format  and  the  0 to  360- 
deg  format  should  be  made  to  aid  with  the  data  interpretation.*  Only  plots 
2 and  3 are  used  for  analysis,  but  plots  1 and  4 can  be  helpful  when  poor  data 
are  encountered.  Interpretation  of  these  plots  can  be  accomplished  with  a 
minimum  background  in  the  analysis  and  trends  of  the  data.  A detailed  example 
of  the  procedures  used  in  interpreting  phase  angle/frequency  plots  is  presented 
in  appendix  B. 

Once  the  plots  have  been  interpreted,  the  phase  angle,  <p,  must  be  determined 
and  recorded  for  certain  frequencies.  (CERF  has  designed  transparent  tem- 
plates which  expedite  this  process.)  The  phase  angle  must  be  continuously 
accumulated  from  0 deg  at  zero  Hz  to  3500  Hz.  On  the  ilSO-deg  plot,  the 
phase  angle  is  read  directly  over  the  frequency  range  for  which  <}>  < 180  deg. 
When  +180  deg  is  reached  (at  the  bottom  of  the  plot)  the  trace  flops  over  to 
the  top  of  the  graph  and  progresses  downward  to  the  bottom  of  the  plot  again. 
The  phase  angle  at  any  frequency  past  the  first  flopovar  is  equal  to  180  deg 

★ ”r 

See  reference  6 for  full  details  on  the  electronic  procedures  for  pro- 
ducing these  graphs. 
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plus  the  angle  represented  by  the  distance  from  the  top  of  the  graph  to  the 
point  on  the  interpretation  curve.  See  figure  20a  for  examples  of  phase 
angles  at  different  points  along  the  curve.  If  a 0 to  360-deg  plot  is  used, 
a phase  angle  of  0 deg  is  at  the  top  of  the  plot  and  360  deg  is  at  the  bottom. 
The  phase  angle  increases  to  360  deg,  at  which  point  a flopover  occurs,  and 
the  trace  begins  again  at  the  top  of  the  plot  (0  deg).  Phase  angles  past  the 
360-deg  flopvoer  point  are  ecf&l  to  360  deg  plus  the  angle  represented  by  the 
distance  from  the  top  of  the  graph  to  the  point  on  the  interpretation  curve. 
See  figure  20b  for  examples  of  phase  angles  at  different  points  along  the 
curve.  The  maximum  phase  angle  depends  on  the  interpretation  of  the  phase 
angle/frequency  plot  and  which  accelerometers  are  used  for  analysis. 


DISPERSION  CURVES 

The  following  equations  are  used  to  compute  the  wavelengths  and  phase  velo- 
cities for  the  dispersion  curve: 


A = -- 


where 


v = fA 

A = wavelength  (ft) 

d = distance  between  compared  accelerometers  (ft) 
i>  = phase  angle  (deg) 
v = phase  velocity  (ft/sec) 
f = frequency  (Hz) 


The  wavelength  equation  was  proven  valid  by  previous  researchers  involved 
with  the  NDPT  van  and  its  development  (ref.  7).  Similarly,  the  equal  ac- 
celerometer spacing  of  2 ft  was  found  adequate  for  most  pavements  tested 
(ref.  26).  A computer  program,  NDTPL0T,  computes  wavelengths  and  phase 


26.  Rao,  H.  A.  B.,  lici vaIUi  of  NontiosLrucLivc  Testa  al  Webb  Air  Faroe  tt/rtc. 

Letter  Report  III,  Civil  Engineering  Research  Facility,  Albuquerque, 
New  Mexico,  August  1 973 - 


velocities  and  plots  the  dispersion  curves.  A sample  problem  is  in- 
cluded in  appendix  B. 

The  interpretation  of  dispersion  curves  so  as  to  select  velocities  for  ma- 
terial property  determination  requires  some  knowledge  of  dispersion  in  a 
layered  media.  A detailed  study  of  dispersion  as  it  relates  to  pavement 
evaluation  is  presented  in  reference  27.  After  the  phase  angle/ frequency 
data  have  been  collected  and  the  dispersion  curve  has  been  plotted,  a 
best  fit  line  is  drawn  through  the  data  points.  Break  points  (abrupt 
changes  in  slope)  on  the  curve  are  used  to  calculate  the  material  proper- 
ties. Pavement  cross-section  data  are  useful  since  there  should  be  a 
branch  of  the  curve  for  each  layer,  including  the  subgrade. 

An  interpretation  of  the  dispersion  curves  for  Test  Pit  6 (Cannon  Air  Force 
Base)  is  shown  in  figure  14  (section  4)  where  a best-fit  curve  through  the 
plotted  points  is  easily  drawn.  The  peak  velocity  occurs  at  5280  ft/sec  and 
represents  the  Rayleigh  wave  velocity  for  the  surface  course.  The  break 
point  of  the  curve  (at  the  end  of  the  high-frequency  data)  occurs  at  1800 
ft/sec  and  represents  the  shear  wave  velocity  of  the  base  course.  It  can 
be  seen  that  the  low-frequency  (long  wavelength)  curve  asymptotically 
approaches  a velocity  of  330  ft/sec  and  represents  the  shear  wave  velocity 
of  the  subgrade  material.  These  velocities  are  then  used  with  either  known 
or  assumed  material  densities  and  Poisson's  ratios  to  determine  the  elastic 
moduli  of  the  material  layers. 


ELASTIC  MODULI 


The  equation  used  to  calculate  the  shear  modulus,  G,  of  the  material  is 


27.  Finn,  Fred,  McCullough,  B.  F.,  Nair,  Keshavan,  and  Hicks,  R.  G., 

/'/  m for  i • r /.  >pm.  uh  of  Non,  l -t'.t  *»/*♦/  .•!<><■  '•/.  •/••:></  f,  <r*  /w.  tw!  i<v.  of 

i-  ■ / of  ' ! ' .iiM’k  <■-.?,  Final  Report  No.  1062- 

2(F),  Materials  Research  and  Development,  Inc.,  Oakland,  California, 
November  1966. 
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where 


V$  = shear  wave  velocity  (ft/sec) 
y = unit  weight  of  material  (lb/ft1) 
g = acceleration  constant  (32.2  ft/sec) 

Young's  modulus,  E,  is  computed  by 

E = 2(1  + y)  G 

where 

v>  = Poisson's  ratio  (usually  assumed  to  be  a reasonable  value) 

The  shear-wave  velocities  for  all  materials  except  the  surface  course  can  be 
inputted  into  these  equations  for  modulus  determinations.  However,  the  Ray- 
leigh wave  velocity  for  the  surface  course  requires  a correction  before  the 
shear-wave  velocity  can  be  computed.  The  equation  used  to  compute  the  shear- 
wave  velocity  from  the  Rayleigh-wave  velocity  is 


where  a is  a function  of  Poisson's  ratio  and  varies  from  0.875  for  v = 0 
to  0.955  for  v = 0.5  (ref.  27). 

When  the  peak  velocities  of  the  dispersion  curve  are  used  to  calculate  elas- 
tic moduli,  very  acceptable  values  are  obtained  for  Portland  cement  concrete; 
however,  high  modulus  values  are  obtained  for  asphaltic  concrete  surfaces. 
Because  of  the  rigid  interparticle  contact  and  bonding  in  Portland  cement  con 
crete,  the  stress  wave  in  such  a layer  propagates  through  the  material  with- 
out any  damping  effects.  However,  in  asphaltic  concrete  pavements,  the 
viscous  damping  properties  of  the  asphalt  binder  at  high  frequencies  retard 
the  phase  angle  shift  and  this  results  in  high  calculated  velocities.  As- 
phaltic concrete  moduli  have  been  observed  by  others  to  increase  with  the 
testing  frequency  (ref.  28).  Therefore,  the  calculated  moduli  must  be  re- 
duced so  as  to  more  accurately  model  the  conditions  in  the  field.  CERF  has 
found  that  a modulus  reduction  factor  of  50  percent  is  reasonable  for  as- 
phaltic concrete  surfaces  of  average  thickness  (3  to  6 in).  A different 


28.  Yoder,  E.  J.,  and  Witczak,  M.  W.,  Principles  oj'  Pavrmnni  Hnniijn,  2nd 
Edition,  John  Wiley  and  Sons,  flew  York,  1975. 


analysis  must  be  performed  when  the  pavement  is  Portland  cement  concrete 
overlaid  with  asphaltic  concrete.  In  this  situation  the  velocity  of  propa- 
gation of  the  Portland  cement  concrete  dominates,  and  the  asphaltic  concrete 
velocity  cannot  be  determined  accurately.  This  is  especially  true  with 
thin  overlays;  here  the  asphaltic  concrete  modulus  must  be  estimated  on  the 
basis  of  the  general  competency  of  the  surfacing. 

Typical  moduli  for  various  types  of  paving  materials,  determined  by  plate- 
bearing tests  performed  at  CEL  (ref.  29),  are  given  in  table  4.  When  moduli 
are  computed  for  the  base  and  subgrade  layers  from  dispersion  curve  data, 
the  values  obtained  are  also  too  high  and  must  be  corrected.  Correction 
procedures  based  on  theoretical  considerations  are  used  to  adjust  the  NDPT 
moduli  to  acceptable  values. 

Two  different  procedures  are  employed  to  correct  the  base  course  moduli 
depending  on  the  type  of  pavement  surface.  A theoretical  approach  is  used 
for  Portland  cement  concrete  pavements.  From  the  theory  of  wave  propagation 
through  layered  media,  the  velocity  should  increase  as  a stress  wave  travels 
from  a lower-speed,  less-dense  material  (i.e.,  base  course)  to  the  higher- 
speed,  more-dense  material  (i.e.,  Portland  cement  concrete  pavement).  The 
governing  equation  for  propagation  of  a shear  wave  through  media  is 


Thus,  the  shear  wave  velocity  in  the  base  course  (layer  2)  is  defined  by 


so  that  the  velocity  increase  as  the  wave  travels  into  the  Portland  cement 
concrete  pavement  (layer  1)  is 


29.  Nielsen,  J.  P.,  !'’nU  i'tnn : Analyfti.n  Hi'lilra'  A 

Technical  Report  R-594,  Civil  Engineering  Laboratory,  Port  Hueneme, 
California,  September  1968. 
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Table  4.  Typical  Moduli  for  Paving  Materials 


Material 

Modulus  Range,  psi 

Portland  Cement  Concrete 

2.0  - 5.0  x 10f> 

Asphaltic  Concrete 

0.1  - 1 .0  x 10* 

Well -Graded  Crushed  Stone  Base 

30  - 90  x 103 

Cement-Treated  Base  (SM,  SC) 

25  - 40  x 103 

Well -Graded  Gravels  (GW) 

Silty  Gravels  (GM) 

10  - 25  x 103 

Silty  Sand  (SM) 

Sandy  Clay  (SC) 

3 - 10  x 103 

Clays  (CL) 

where 

V'  = shear  velocity  in  the  base  course  as  indicated  by  the  dispersion 
curve 

Vs  = actual  in-situ  shear  velocity  in  the  base  course 

G,  = shear  modulus  of  the  surface  layer 

G'  = shear  modulus  of  the  base  course  as  derived  from  V' 

s? 

The  above  equation  can  be  solved  for  V<-  , and  thus  the  in-situ  shear  modulus 

5 2 

can  be  determined  from  8 = V<-  p.  This  analysis  was  performed  on  the  Cannon 
Air  Force  Base  data  and  yielded  a corrected  velocity  range  of  625  to  325  ft/ 
sec  for  a stabilized  base  course.  These  velocities  correspond  to  Young's 
moduli  of  30,000  to  50,000  psi,  which  are  satisfactory  values  for  a base 
course. 
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The  following  sample  calculation  is  performed  with  data  from  Test  Pit  2 at 
Cannon  Air  Force  Base: 


VR  = 6200  ft/ sec  (surface  Rayleigh  wave,  a = 0.90) 


V$  = 6889  ft/ sec 

i 

y = 145  lb/ft3 

o = 0.20 

i 


R,  “ f6889^  (^32.2 J 144 J 1,48  x 1Q6  psi 
G;  = ( 21 30 ) 2 (^5-4)  = 1-32  x 10s  psi 


= 2130  ft/sec  (disper- 
2 sion  curve  data) 

Y = 135  lb/ft3 

2 

v = 0.30 
2 


1 (1.32  x 105] 

1 1 

;i45) 

(Tl.48  x 10O 

1 1 

T35l 

<%  - <669>2  (srr)  (m) 


¥.7)  hir) ' ,2-65“  Psi 


E = 32,900  psi 

? 


An  alternate  and  simpler  approach  is  to  use  the  approximate  equation 


(E')2 


where 


E = adjusted  Young's  modulus 

E^  = Young's  modulus  obtained  from  the  dispersion  curve  shear-wave 
velocity  for  the  base  course 
E^  = Young's  modulus  of  Portland  cement  concrete 

When  this  approach  is  used,  the  following  results  are  obtained: 


EJ  = (1.32  x 105)  (2)  (1  + 0.30)  = 3.43  x 105  psi 


E = (1.48  x 106)  (2)(1  + 0.20)  = 3.55  x 10r'  psi 
1 


(3.43  x 10T  „ 

13. 55  T T0TT  = 33  ’ 160  PS1 


IUT,'"U,U'- 


" ~ ' 


u"'"nut  11 


When  the  pavement  surface  is  asphaltic  concrete,  a 50-percent  reduction 
is  applied  to  the  base  course  modulus  value.  Although  past  experience  with 
asphaltic  concrete  pavements  has  indicated  that  this  factor  is  appropriate, 
no  theoretical  development  can  be  presented  to  validate  this.  However,  as 
with  the  correction  procedure  for  the  surface  modulus,  it  appears  that  the 
viscous  damping  properties  of  the  asphaltic  concrete  interfere  with  the  de- 
termination of  the  actual  velocities  of  propagation.  An  empirical  correction 
factor,  therefore,  is  deemed  necessary. 

Subgrade  modulus  values  are  corrected  based  on  an  empirical  correction  factor 
of  50  percent  for  both  asphaltic  concrete  and  Portland  cement  concrete  pave- 
ments. This  value  was  arrived  at  through  the  studies  at  Cannon  and  Laughlin 
Air  Force  Bases.  When  applied  to  the  data,  the  subgrade  moduli  were  corrected 
to  an  acceptable  range  of  3,000  to  10,000  psi.  A summary  of  the  modulus  cor- 
rection procedure  is  presented  in  table  5. 


Table  5.  Modulus  Correction  Procedure 


Pavement  Layer 

Correction 

Surface 

Portland  Cement  Concrete 
Asphaltic  Concrete 

Use  actual  value  determined. 
Reduce  by  50  percent. 

Basp/ Subbase  with: 
PCC  Surface 
AC  Surface 
AC/PCC  Surface 

Corrected  E = E ?/E 

2 2 i 

Reduce  by  50  percent 
Reduce  by  50  percent 

Subgrade  with: 
PCC  Surface 
AC  Surface 
AC/PCC  Surface 

Reduce  by  50  percent 
Reduce  by  50  percent 
Reduce  by  50  percent 
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APPENDIX  A 


OPERATIONAL  INSTRUCTIONS 
FOR 

LABORATORY  DATA-REDUCTION  INSTRUMENTATION 

INTRODUCTION 

This  appendix  provides  brief  instructions  for  the  operation  of  the  laboratory- 
based,  data-reduction  instrumentation  shown  in  figure  A1 . This  equipment  is 
similar  to  that  installed  in  the  NDPT  van.  These  instructions  have  been 
written  with  the  understanding  that  operators  of  this  equipment  are  thoroughly 
familiar  with  the  equipment  in  the  van;  hence,  these  instructions  lack  the 
details  provided  in  the  user's  manual  for  the  van  (ref.  6). 


INITIAL  SETUP 


(1)  Turn  on  the  main  power  (lower  left  section  of  rack).  Indicator 
lights  above  the  left  and  right  racks  should  be  lit. 

(2)  Clean  the  tape  deck  record/reproduce  heads  with  cleaning  fluid  accord- 
ing to  the  instructions  and  recommendations  included  in  the  manufac- 
turer's operation  manual.  Clean  heads  with  Ampex  head  cleaner.  Clean 
capstan  (puck)  with  denatured  alcohol  only. 

(3)  Turn  on  the  tape  deck  power. 

(4)  Install  a scratch  tape  in  the  tape  deck  and  depress  the  STOP  READY 
Button. 

(5)  Insure  that  other  equipment  is  receiving  power  and  that  power  switches 
are  on. 

(6)  Verify  that  the  power  supply  for  the  punch  panel  indicator  lights 
(lower  right  section  of  rack)  registers  12  Vdc.  Adjust  if  necessary. 

Allow  a 30-min  warmup  period  before  performing  any  calibrations.  During 

this  time,  field  notes  and  audio  tape  documentations  should  be  reviewed  to 

determine  the  exact  channels  to  be  used  to  record  the  field  data. 


EQUIPMENT  CALIBRATION 
Differential  Amplifiers 

All  necessary  calibration  procedures  are  outlined  in  detail  in  the  manufac- 
turer's manuals.  However,  once  these  calibrations  have  been  completed,  the 
following  must  be  performed  before  other  system  calibrations  are  made: 

(1)  Turn  the  AMPLIFIER  OUTPUT  Switches  (on  punch  panel)  for  Channels  1 
and  2 to  FILTERED. 
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(2)  Turn  the  REFERENCE  FREQ,  and  ANALOG  FREQ.  Switches  to  DIRECT. 

(3)  Turn  the  FILTER  Switches  for  both  amplifiers  to  10  kHz. 

(4)  Turn  the  VAR  GAIN  (variable  qain)  Switches  for  both  amplifiers  counter- 
clockwise to  DETENT  (calibrate). 

Phase  Computer 

(1)  Remove  softwire  connections  to  the  phase  computer  inputs  for  Channels 

A and  B before  calibration. 

(2)  Push  the  METER  RANGE  Switch  to  + 180  . 

(3)  Push  the  FILTER  SELECTION  Switch  to  1 (latched).  This  indicates  a 
minimum  amount  of  signal  smoothing. 

(4)  Turn  the  N0RMAL/180  Switch  to  NORMAL  (released). 

(5)  Turn  the  3mV-5V/xlOO  Switch  to  3mV-5V  (released)  for  both  channels. 

(6)  With  no  signal  applied,  and  adjusting  the  PHASE  ZERO  Control  Potentio- 

meter, correct  the  digital  phase  indicator  to  read  + 180.0". 

(7)  Depress  the  0-360  METER  RANGE  Switch  and  verify  that  the  indicator 
reads  360.0°. 

NOTE 

An  exact  reading  of  360'  may  not  show  on  the  indicator; 
the  reading  may  deviate  a few  tenths  of  a degree  depend- 
ing on  the  + 180"  calibration. 

(8)  Return  the  METER  RANGE  Switch  to  + 180". 

(9)  Replace  softwire  connections  to  the  phase  computer  inputs  for  Channels 
A and  B. 


J 
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Carrier  Generator 


NOTE 

Refer  to  the  manufacturer's  manual  for  laboratory  calibration. 
(1)  Turn  the  OPER  CAL  Switch  to  CAL. 

(?)  Verify  that  the  needle  is  in  the  green  area  on  the  tuning  amplitude 
meter.  If  it  is  not,  refer  to  the  manufacturer's  manual  for  the  pro- 
per calibration  procedure. 


NOTE 

Periodically  during  the  test,  check  that  the  tuning  am- 
plitude meter  needle  is  reading  in  the  green  area;  adjust 
if  necessary  using  the  TUNING  AMPL  (tuning  amplitude) 

Potentiometer  with  the  OPER  CAL  Switch  on  OPER. 

Tracking  Filter 

(1)  Verify  that  the  OPER  CAL  Switch  on  the  carrier  generator  is  on  CAL. 

(?)  Turn  the  RANGE  Switch  for  Channels  1 and  2 to  CAL. 

(3)  Turn  the  RELATIVE  RANGE  Vernier  for  both  channels  to  CAL. 

(4)  Turn  the  METER  MULTIPLIER  Switch  for  both  channels  to  CAL. 

(5)  Verify  that  the  FILTER  SELECTION  Switch  on  the  phase  computer  is  on  1. 

(6)  Turn  the  MODE  SELECTION  Switch  to  SINE  AVG  (SEC)  1. 

(7)  Using  the  FREQ  Trim  Potentiometer  of  the  carrier  generator,  adjust  for 
maximum  signal  strength  of  both  channels  (indicated  by  the  signal 
strength  meters  of  the  tracking  filter). 
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NOTE 


It  may  be  necessary  to  adjust  one  of  the  channels  to  a 
maximum  signal  and  fine  tune  the  other  as  indicated  in 
step  8. 

(8)  Using  the  CAL  Trim  Potentiometer  for  Channel  1 Filter  1,  adjust  the 
Channel  1 signal  meter  to  read  10  (full  range). 

(9)  Repeat  step  8 using  Channel  2 Filter  1 CAL  Trim  Potentiometer  and 
Channel  2 signal  meter. 


(10)  Turn  the  carrier  generator  OPER  CAL  Switch  to  OPER. 


(11)  Turn  the  RANGE  Switch  for  Channels  1 and  2 to  3Vrms. 

(12)  Turn  the  METER  MULTIPLIER  Switches  of  the  tracking  filter  to  1. 

(13)  Input  a calibration  signal  to  both  channels  to  adjust  the  tracking 

filter.  The  Spectral  Dynamics  Sweep  Oscillator  Servo  (upper  left 

rack)  is  used  for  this  purpose.  The  SOS  signal  is  hardwired  to  the 

patch  panel  CAL  SIGNAL  OUT  Terminal  and  softwire  connections  are  made 
to  input  this  signal  to  the  tracking  filter  (fig.  A2).  Adjust  the 
sweep  oscillator  as  follows: 

(a)  Set  FUNCTION  Switch  to  SINE  WAVE. 

(b)  Set  LEVEL  Potentiometer  to  6 (relative  gain  setting). 


(c)  Switch  MODE  SELECTION  Switch  to  LINEAR  SWEEP  SINGLE. 


(d)  Switch  Multiplier  A to  10. 

(e)  Set  Multiplier  B to  10. 

(f)  Set  sweep  rate  (30  to  50  Hz/sec  is  recommended) . 


(g)  Sweep  up  to  3500  Hz  and  hold  at  this  frequency. 
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(14)  Adjust  the  TUNING  AMPL  Trim  Potentiometer  of  the  carrier  generator 
until  the  meter  reads  mid-scale  in  the  green  area. 

(15)  Set  the  GAIN  Switches  of  the  differential  amplifiers  so  that  the  signal- 
strength  meters  of  the  tracking  filter  indicate  approximately  the  same 
mid-scale  reading. 

NOTE 

The  GAIN  Switches  should  be  set  at  identical  positions. 

(16)  Adjust  the  Channel  1 Filter  1 Trim  Potentiometer  on  the  tracking  filter 
to  indicate  a system  phase  angle  of  0"  (+-  0.2  ) on  the  phase  meter. 

AMPEX  PR-2200  Tape  Recorder/Reproducer 

For  calibration  procedures  refer  to  the  manufacturer's  operation  manual. 
Knowledge  of  the  quality  level  of  recording  is  necessary  to  properly  cali- 
brate and  adjust  the  reproduce  amplifiers  for  data  reduction.  Only  a 
qualified  technician  should  perform  the  calibrations  given  in  the  manufac- 
turer ' s manual . 

The  tape  deck  is  equipped  with  special  reproduce  heads  which  can  be  adjusted 
to  correct  for  skew,  which  could  introduce  errors  into  the  phase  angle  com- 
putation. The  following  procedure  is  required  to  adjust  the  t^pe  deck  skew 
between  different  channels  for  accurate  phase-angle  determinations  of  field 
data . 

(1)  Press  the  TAPE  SPEED  Switch  (i.e.,  7 1/2  in/sec). 

(2)  Disconnect  CAL  SIGNAL  OUT  from  amplifier  input  terminals  (AMP  IN  1 and 
AMP  IN  2)  on  patch  panel  and  connect  to  tape  deck  channel  inputs  which 
correspond  to  Accelerometers  1 and  5 of  the  field  data  as  shown  in 
figure  A3. 
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Figure  A3.  Supplement  Wiring  for  Calibration  of  Tape  Deck  Skew 
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(3)  Connect  Channels  1 and  9 tape  output  to  amplifier  input  AMP  IN  1 and 
AMP  IN  2,  respectively. 

(4)  Turn  the  REFERENCE  FREQ  and  ANALOG  FREQ  Switches  on  the  punch  panel  to 
PLAYBACK. 

(5)  Place  tape  deck  in  record  forward  mode  by  simultaneously  depressing  the 
RECORD  and  FWD  Buttons. 

(6)  Adjust  the  TUNING  AMPL  Trim  Potentiometer  on  the  carrier  generator  until 
the  TUNING  AMPLITUDE  Meter  is  in  the  green  area. 

CAUTION 

Adjustments  should  be  very  minor;  do  not  overextend  the 
range  of  the  skew  adjusting  screw.  Damage  to  the  tape 
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heads  can  result  if  the  screws  are  extended  past  their 
normal  range. 


(7)  While  monitoring  the  digital  phase  indicator,  adjust  the  skew  adjusting 
screw  for  the  odd-numbered  tape  channels  to  minimize  the  skew  error  be- 
tween channels  indicated  by  the  phase  reading. 

NOTE 

The  adjusting  screws  are  labeled  0 and  E for  the  odd-and 
even-numbered  tape  heads,  respectively.  The  digital  phase 
indicator  should  show  a phase  angle  of  0°  when  calibrated 
properly. 

(8)  Check  other  channel  combinations  on  the  same  head  to  verify  calibrations 

NOTE 

The  same  procedure  must  be  performed  for  the  even- 
numbered  channels  with  the  E skew  adjusting  screw. 

(9)  Stop  the  tape  deck. 

(10)  Disconnect  the  following: 

(a)  Signal  generator  output  from  tape  deck  inputs  (i.e..  Channels  1 
and  9)  on  the  patch  panel. 

(b)  Signal  generator  output  from  carrier  generator  (Channel  10)  on 
the  patch  panel . 

(c)  Generator  signal  from  CAL  SIGNAL  OUT  on  the  patch  panel. 

(11)  Turn  off  signal  generator. 


(12)  Remove  scratch  tape. 


(13)  Install  field  test  data  tape. 


X-Y-Y'  Plotter 

The  X-Y-Y'  plotter  is  capable  of  recording  two  different  signals  on  the  Y- 
axis.  This  capability  is  used  to  calibrate  the  plotter  and  to  create  two 
different  plots.  One  axis  is  used  to  plot  the  + 180"  plot  and  the  other  is 
used  to  plot  the  0-360'  plot.  By  using  the  two  different  axis  of  the  plotter, 
recalibration  of  the  plotter  is  unnecessary  when  both  plots  are  desired. 

(1)  Place  plot  paper  on  plotter  and  turn  the  CHART  Switch  to  HOLD. 

(2)  Properly  install  plotter  pens  in  pen  carriages  (see  instruction  manual) 
and  turn  the  PEN  Switch  to  LIFT. 

(3)  Turn  RESPONSE  Switches  (3)  to  SLOW. 

(4)  Turn  on  SERVO  Switch. 

Y-l  Axis 

The  Y-l  axis  of  the  plotter  is  used  to  create  the  + 1 80c  plot.  The  follow 
ing  calibration  of  the  Y-l  axis  is  required  for  this  plot. 

(1)  With  no  voltage  applied  to  the  Y-l  input,  adjust  the  pen  to  regis- 
ter on  the  center  of  the  plot  using  the  ZERO  Adjust  knob. 

(2)  Observing  standard  signal  polarity,  input  + 1.8  Vdc  from  a preci- 
sion power  supply  (e.g.,  DIAL-A-VOLT)  to  the  Y-l  axis  of  the 
plotter. 

(3)  Turn  the  Y-l  RANGE  SELECTOR  Switch  to  .lV/in. 

(4)  Turn  the  POLARITY  Switch  to  -UP. 

(5)  Using  the  Y-l  RANGE  CAL  Vernier,  adjust  the  Y-l  axis  to  read  + 180" 
on  the  plot  (i.e.,  bottom  of  the  plot). 


(6)  Flip  the  POLARITY  Switch  to  HIP  and  verify  that  the  plotter  in- 
dicates - 180°  (i.e.,  top  of  the  plot). 

(7)  After  calibration,  turn  the  POLARITY  Switch  to  -UP.  (This  will 
result  in  a plot  of  + 180  at  the  bottom  of  the  plot  and  - 180'J 
at  the  top. ) 

(8)  Disconnect  calibration  input  voltage. 

Y-2  Axis 

The  Y-2  axis  of  the  plotter  is  used  to  create  the  0-360'J  plot.  The 

following  calibration  of  the  Y-2  axis  is  required  for  this  plot. 

(1)  With  no  voltage  applied  to  the  Y-2  input,  adjust  the  axis  using  the 
Y-2  ZERO  Adjust  Knob  of  the  plotter  to  correspond  to  0 at  the  top 
of  the  plot. 

(2)  Observing  standard  signal  polarity,  input  + 3.6  Vdc  from  a preci- 
sion power  supply  to  the  Y-2  axis  of  the  plotter. 

(3)  Turn  the  Y-2  RANGE  SELECTOR  Switch  to  . lV/in. 

(4)  Turn  the  Y-2  POLARITY  Switch  to  -UP. 

(5)  Using  the  Y-2  RANGE  CAL  Vernier,  adjust  the  plot  to  360  (i.e.,  bot 

tom  of  plot) . 

(6)  Disconnect  calibration  input  voltage. 


(7)  Reconnect  phase  output  ,o  appropriate  axis. 


X-Axis 


(1)  Set  the  RANGE  SELECTOR  Switch  of  the  X-axis  to  an  appropriate  range 
for  the  recorded  test  data. 

(2)  Turn  the  POLARITY  Switch  to  +RT  (right).  (This  will  plot  with  the 
frequency  increasing  to  the  right.) 

(3)  Proceed  to  low  end  of  test  (i.e.  , 10  Hz). 

(4)  Using  the  analog  frequency  from  the  data  tapes,  set  the  lower  fre- 
quency limit  (i.e.,  10  Hz)  with  the  ZERO  Adjust  Knob  of  the  X-axis. 

(5)  Proceed  to  the  beginning  of  the  test  (i.e.,  350  or  3500  Hz). 

(c*)  Using  the  RANGE  CAL  Vernier  of  the  X-axis,  adjust  the  upper  fre- 
quency limit  (i.e.,  350  or  3500  Hz). 

(7)  Turn  the  three  RESPONSE  Switches  to  FAST. 

(8)  Rewind  tape  to  beginning  of  recorded  field  test  data. 

PLOTTING  OF  FIELD  DATA 

(1)  Connect  phase  angle  output  on  patch  panel  to  Y-l  input  of  plotter  for  a 

* 180  plot.  (If  a 0-360'’  plot  is  desired,  the  phase  angle  output  should 
be  connected  to  the  Y-2  input.  ) 

(2)  Verify  that  the  MET'R  RANGE  Switch  of  the  phase  computer  is  on  + 180’. 

(If  a 0-3GU'  plot  is  desired,  the  METER  RANGE  Switch  shoul.i  be  on 
0-360  . ) 

(3)  Verify  that  the  N0RMAL/180  Switch  is  on  NORMAL  (released). 

(4)  Verify  that  the  3mV-5V/xlOO  Switch  is  on  3mV-5V  (released). 
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(5)  Verify  that  the  X-axis  is  connected  to  the  analog  frequency  tape  output 
channel  (i.e.,  Channel  4). 

(6)  Turn  IRIG  MODE  SELECTOR  Switch  to  TRANS. 

I 

(7)  Turn  SYNC  Switch  of  IRIG  generator/reader  to  IRIG  B. 

(8)  Turn  INPUT  Switches  on  the  IRIG  qenerator/reader  to  + and  FWD. 

NOTE 

1 he  position  of  the  FILTER  Switch  of  the  IRIG  generator/ 
reader  is  determined  by  the  quality  of  the  recorded  IRIG 
signal . 

(9)  Connect  accessory  equipment  for  data  monitoring  (i.e.,  frequency  coun- 
ter, oscilloscope,  etc.)  to  proper  channels,  utilizing  the  tape  out- 
puts . 


(10)  Verify  that  the  following  switches  located  on  the  punch  panel  are  in  the 
indicated  positions: 


Swi tch 


Posi tion 


Panel  Power 
Analog  Frequency 
Reference  Frequency 
Amplitude  1 Output 
Amplitude  2 Output 
IRIG  Reader 


PLAYBACK 

PLAYBACK 

FILTERED 

FILTERED 


(11)  Connect  desired  accelerometer  channels  to  proper  amplifier  input. 

(The  inboard  accelerometer  should  be  connected  to  the  AMP  IN  1 Terminal; 
the  outboard  accelerometer  should  be  connected  to  the  AMP  IN  2 Terminal.) 

(12)  Insure  that  the  SERVO  Switch  on  the  plotter  is  on  STAND  BY. 


(13)  Play  back  data  in  the  forward  inode  by  pressing  the  FWD  Button  on  the 
recorder. 


(14)  Check  tracking  filter  signal  level  and  adjust  gain  levels  to  pro- 
vide mid-scale  readings. 

(15)  Flip  the  SERVO  Switch  on  the  plotter  to  ON. 

(16)  Turn  the  PEN  Switch  to  RECORD  when  frequency  sweep  begins. 

(17)  When  the  frequency  sweep  is  completed,  turn  PEN  Switch  to  LIFT. 


(18)  Flip  the  SERVO  Switch  to  STAND  BY. 

(19)  If  other  accelerometer  data  are  desired,  rewind  the  tape  to  beginning 
of  test  and  connect  the  desired  accelerometer  channels  to  the  AMP  IN 
Terminals  and  repeat  steps  2,  and  11  through  18.  If  data  from  another 
test  are  desired,  proceed  to  next  test  tape  position  and  repeat  steps 
2,  and  11  through  18. 


SHUTDOWN 


(1)  Rewind  the  data  tape  and  remove  it  from  the  tape  deck. 

(2)  Turn  off  the  tape  deck  power  and  allow  the  internal  timing  sequence 
to  turn  off  the  tape  deck. 


(3)  Turn  off  instrumentation  rack  power  and  accessory  equipment. 
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APPENDIX  B 

SAMPLE  DATA-REDUCTION  PROBLEM 

This  appendix  outlines  the  recommended  procedure  for  reducing  the  phase 
angle/frequency  plots,  constructing  the  dispersion  curves,  and  determining 
modulus  values.  Actual  field  data  from  Shaw  Air  Force  Base  are  utilized 
in  the  example. 

The  actual  phase  angle/frequency  plots  for  Shaw  are  shown  in  figures  B1  and 
B2  for  the  low-frequency  and  high-frequency  sweeps,  respectively.  The  gen- 
eral trend  ( interpretation)  of  the  plots  is  shown  by  the  heavy,  solid,  smooth 
line.  In  figure  B1 , the  trend  of  the  data  can  easily  be  seen  and  no  problem 
exists  in  interpretation  above  125  Hz.  However,  below  125  Hz,  the  data  show 
no  trend  and  the  interpretation  shown  was  arrived  at  by  examination  of  other 
data  at  Shaw  and  through  previous  experience  with  phase  angle/frequency  data. 
The  trend  from  point  A through  point  E is  the  appropriate  interoretation.  An 
interpretation  based  on  a straight  line  from  the  zero  position  to  180  deg 
at  125  Hz  is  not  accurate  and  contrary  to  trends  observed  from  other  NDPT 
data.  It  should  be  pointed  out  that  the  interpretations  of  the  data  must  be 
j'opoad.  through  a phase  angle  of  180  deg  or  multiples  thereof.  The  actual 
field  data  never  plot  to  exactly  130  deg.  This  is  due  to  the  effect  of  the 
sampling  rate  of  the  phase  computer.  The  sampling  rate  could  be  increased 
(i.e.,  make  shorter  time  intervals  between  phase  computations)  but  the 
amount  of  noise  (illustrated  by  the  erratic  traces)  would  also  be  increased. 
Therefore,  the  data  were  forced  through  180  deg  at  125  Hz.  Below  this  fre- 
quency, the  interpretation  represents  a typical  pavement  response.  Figure  B1 
may  be  used  as  a model  for  the  general  trend  of  the  low-frequency  data,  in- 
dependent of  pavement  type.  Specific  points  on  the  interpretation  curves 
are  labeled  with  the  corresponding  phase  angles. 

The  phase  angle/frequency  coordinates  of  points  A through  E can  be  read 
directly  from  the  plot  (fig.  Bl);  those  of  points  F through  I can  also  be  read 
directly  from  the  plot  (fig.  B2);  but  points  J through  ? require  some  expla- 
nation. If  the  interpretation  line  is  followed  beginning  at  0 Hz  (fig.  B2), 
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Accelerometers  1 and  5 


Figure  Bl.  Low-Frequency  Phase  Angle/Frequency  Data  for  Shaw  Air  Force  Base 


Figure  B2.  High-Frequency  Phase  Angle/Frequency  Data  for  Shaw  Air  Force  Base 


the  phase  angle  increases  to  180  deg  at  a frequency  of  912  Hz,  where  a flopover 
occurs.  The  phase  angle  then  accumulates  continuously.  For  example,  point  M 
is  180  deg  plus  the  phase  angle  represented  by  the  distance  M-M'  (i.e.,  19  deg). 
The  total  phase  angle  at  point  M is  199  deg.  Similarly,  point  N and  point  0 
represent  total  phase  angles  of  410  and  483  deg,  respectively.  Another  flop- 
over  occurs  at  a frequency  of  2687  Hz;  the  corresponding  phase  angle  is  540 
deg.  The  phase  angle  for  point  P is  determined  by  adding  the  angle  P-P’  (i.e., 
140  deg)  to  540  deg,  or  680  deg.  CERF  has  developed  transparent  overlays  which 
can  be  placed  over  the  plots  to  expedite  the  data  reduction.  The  coordinates 
of  the  interpreted  plot,  that  is  the  phase  angle/frequency  values,  are  taken 
every  25  Hz  from  25  to  125  Hz  for  the  low-frequency  sweep  and  every  125  Hz  from 
250  to  3500  Hz  for  the  high-frequency  sweep.  Note  carefully  that  accelerometers 
1 and  5 or  2 and  5 are  used  for  the  former;  accelerometers  1 and  3 or  2 and  3 
are  used  for  the  latter.  It  is  a simple  matter  to  prepare  a work  sheet  onto 
which  these  data  can  be  recorded. 

When  these  data  have  been  collected  from  the  low-and  high-frequency  sweep 
plots,  the  computer  program  NDTPL0T  (appendix  C)  or  manual  computation  of 
wavelength  and  velocity  is  required  to  develop  the  dispersion  curve.  The 
equations  are  as  follows: 

. 360d 

A - 

V = fX 


where 

X = wavelength  (ft) 

d = spacing  between  accelerometers  (ft) 

<j>  = phase  angle  (deg) 
v = phase  velocity  (ft/ sec) 
f = frequency  (Hz) 

It  should  be  pointed  out  that  the  wavelength  equation  contains  the  parameter 
d,  accelerometer  spacing,  and  that  the  phase  angle/frequency  data  were  obtained 
from  two  plots  and  two  different  accelerometer  spacings.  If  the  manual  compu- 
tation method  is  used,  the  appropriate  values  of  d should  be  used.  This  means 
that  d = 2 ft  for  the  high-frequency  (250  to  3500  Hz)  data  and  d = 4 ft  for  the 
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low-frequency  (25  to  125  Hz)  data.  However,  with  NDTPLOT,  the  phase  angle  for 
the  low-frequency  data  is  divided  by  two  and  then  inputted  to  the  program  along 
with  the  high-frequency  data.  A spacing  of  2 ft  is  then  inputted  to  the  pro- 
gram because  the  program  is  incapable  of  using  different  accelerometer  spacings 
for  two  sweeps  placed  on  the  same  plot.  Thus,  it  is  necessary  to  divide  the 
actual  phase  angle  by  two  for  the  low-frequency  data  only. 

The  phase  angle/ frequency  data  and  the  resulting  dispersion  curve  are  shown 
in  figures  B3  and  B4,  respectively.  Figure  B3  illustrates  the  format  used  to 
input  the  phase  angles  for  selected  frequencies  to  the  NDTPLOT  program.  The 
frequencies  are  inputted  to  the  program  by  a data  statement  incorporated  in 
the  program.  Therefore,  care  must  be  taken  to  insure  that  the  phase  angle 
values  are  positioned  on  the  Fortran  coding  form  so  that  they  correspond  to 
the  proper  frequencies  (appendix  C).  Notice  that  the  phase  angles  indicated 
in  figure  B1  corresponding  to  frequencies  of  25,  50,  75,  100,  and  125  Hz  have 
been  divided  by  two  on  the  coding  form. 


IBM  FORTRAN  Coding  Form  o.,  .•  ja 

[ - ' w>T  pL ot  pr  j ~ ~j~ 
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Figure  B3.  Sample  NDTPLOT  Fortran  Coding  Form 
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5750  ft/sec 


Figure  B4.  Composite  Dispersion  Curves  for  Shaw  Air  Force  Base 


Interpretation  of  the  dispersion  curve  can  easily  be  made  and  the  critical 
velocities  can  be  obtained.  The  Rayleigh  peak  velocity  of  the  surface 
course  of  the  Portland  cement  concrete  pavement  is  5750  ft/sec;  the  shear 
wave  velocity  of  the  base  course  is  1550  ft/sec;  and  the  subgrade  shear 
wave  velocity  is  360  ft/sec.  Next,  the  uncorrected  moduli  are  calculated 
by 


where 

G = shear  modulus  (lb/ft2) 

V$  = shear  wave  velocity  (ft/sec) 

Y = unit  weight  (lb/ft3) 
g = acceleration  constant  (32.2  ft/sec/sec) 

The  Rayleigh  velocity  of  5750  ft/sec  must  be  converted  to  a shear  velocity 
by  dividing  by  0.9  (used  for  both  Portland  cement  concrete  and  asphaltic  con 
Crete).  This  yields  a shear  wave  velocity  of  6389  ft/sec,  which  corresponds 
to  a Young's  modulus  of  3.06  x 106  psi.  A summary  of  the  calculations  is  as 
follows : 

Uncorrected 
Modul i 

Vs>  ft/sec  y,  lb/ft3  _v G,  psi  E,  psi 


Portland  Ce 


ment  Concrete 

6389 

145 

0.25 

1.28 

X 

106 

3.19 

X 

106 

Base  Course 

1550 

120 

0.35 

6.22 

X 

104 

1.68 

X 

105 

Subgrade 

360 

105 

0.38 

2.93 

X 

103 

8.10 

X 

103 

The  modulus  for  the  base  course  is  very  high  when  the  dispersion  curve  veloc 
ity  of  1550  ft/sec  is  used.  However,  the  modulus  can  be  corrected  with  the 
correction  procedure  presented  in  section  5. 


E = 
2 


where 

E2  = corrected  Young's  modulus  (psi) 
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E^  = uncorrected  Young's  modulus  based  on  dispersion  curve  velocity  (psi) 
E = Young's  modulus  of  surface  course  (psi) 


E 

2 


(1.68  x 105)? 
13719  x 10*) 


= 8847  psi 


The  subgrade  modulus  is  reduced  50  percent  from  the  calculated  value  (i.e., 
8100  psi)  to  a corrected  value  of  4050  psi. 


The  structural  analysis  of  the  pavements  is  accomplished  by  implementing 
PREDICT.  Appendix  D contains  a user’s  manual  for  PREDICT. 


APPENDIX  C 

PROGRAM  NDTPLOT  USER'S  MANUAL 


NDTPLOT  USERS  MANUAL 
CERF  VERSION  15  APRIL  1976 
WRITTEN  FOR  AFCEC  REDUCTION  OF  NDPT  DATA 


A CARD  FORMAT  (15) 

COLUMN  VARIABLE  DESCRIPTION 
1-5  NOTS  TOTAL  NUMBER  OF  TEST  STATIONS 
B CARD  FORMAT  (4A10,3X,I2) 


1 


i 


I 


COLUMN 

VARIABLE 

DESCRIPTION 

1-20 

TITLE 

NAME  OF  FACILITY  (EG. NELSON  AFB) 

21-40 

RUNWAY  OR  TAXIWAY  NUMBER  AND  STATION  OR  TEST  POINT 

44-45 

I PAIRS 

NUMBER  (EG. RW  03R-21L  STA. 41+00) 

NUMBER  Or  ACCELEROMETER  PAIRS  PER  TEST  STATION 

(COL. 45  = 1) 

C CARD  FORMAT  (12 ,3X,2A10,15 ,F5.0) 

COLUMN 

VARIABLE 

DESCRIPTION 

1-2 

IAPAIR 

ACCELEROMETER  PAIR  DESIGNATION  NUMBER 

6-25 

TITLEI 

(COL. 2 - 1) 

CAN  BE  USED  TO  DOCUMENT  TEST  DATA 

26-30 

NMPT 

(WILL  APPEAR  ON  PLOT  TITLE) 
NUMBER  OF  PHASE  ANGLE  POINTS 

31-35 

SPACE 

(MUST  CORRESPOND  TO  THE  NUMBER  OF  POINTS  INPUT  ON 
(D  CARDS.  COL. 29-30  = 32) 

SPACING  BETWEEN  ACCELEROMETERS 

(COL. 35  = 2) 

D CARD  FORMAT  (8F10.0) 

COLUMN 

VARIABLE 

DESCRIPTION 

1-10 

PHI 

PHASE  ANGLE 

11-20 

21-30 

ETC. 

(THE  PHASE  ANGLES  INPUT  ON  THE  D CARDS  MUST  CORRESPOND  TO  THE 

(FOLLOWING  FREQUENCIES  25,50,75,100,125,250,375,500,625,750, 
(875,1000,1125,1250,1375,1500,1625,1750,1875 ,2000,2125,2250, 
(2375 ,2500 ,2625 ,2750,2875 ,3000 ,3125 ,3250 ,3375 ,3500. 

(IF  ANY  PHASE  ANGLE  IS  OMITTED  FOR  ANY  OF  THE  ABOVE  FREQUENCIES 
(A  PHASE  ANGLE  =0.  MUST  BE  ENTERED  FOR  THOSE  FREQUENCIES  IN  THE 
APPROPRIATE  COLUMNS.  FOUR  D CARDS  ARE  REQUIRED  TO  INPUT  ALL 
(PHASE  ANGLES  FOR  THE  ABOVE  FREQUENCIES.  IE.  32  VALUES 
(B,C,D  CARDS  MUST  BE  REPEATED  FOR  NOTS  NUMBER  OF  TEST  STATIONS 
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p 


program  NOTPLOT  

***********  o^?PM  Fnw  THF  are  FORCE  WEAPONS  LABORATORY  * 

« THIS  PROGRAM.  wRITTFN  FOR  I^thFRFO  WITH  THF  NON-OF.STRIJCT  I VF  * 

*RY  AMN.  P.  H.  ^AKFR.  USES  DATA  GATHER  O WM  VELOCITY* 

•PAVEMENT  TEST  APPARATUS  TO  CALJULATF  THE  WAVM SECTION,  THF* 

* OF  A SHOCK  w*y^TJH^^INC\\,FcU0N  RECTILINEAR  PAPER  TO  OBTAIN  A* 
^tsSfrsio^cuSvPI \J!f  prints1'^0  Vues  of  wavelength  ano  velo  : 

!2rr  ...  . V.'mil  to  CAUlJ 


.&ly'i?RMS5"~B»E3?3lD  AY  RILL  “On-E  R.Y  1975.  J 

bn : si  msmMpB 

VjWcU  l K8E-SRFtSKl«?  Sppfa^on^Jot^titlf  ; 

KSSJf  ; gy2c?S«%E?SJ£NP?i5T^8f?C«0-FTt-W«i  USED  I 

*FRFO  = FREQUENCY  OF  VTHRATIOM  • 

•PHI  = PHASE  ANGLE  • 

•h AVI GT  = WAVFLENGTH  • 

:s5^*;w********2*in*******************r**t*(*:tc:r* 

n T MF WG  T ON  FRFQ ( 1 SO ) , Ph  I (ISO)  «WAVLGT  (ISO)  » W A VVF  L ( I SO ) 

DIMENSION  MAXWLH?)  • NAXWVM2)  !?II?WL,1?)  *mINWV(1?) 
n I MENS I ON  LGT1 (6« ISO)  .VE L 1 ( 6. ISO ) 

RISESISS  ?rTtI{i??l?ftsiTif:SU«*. 

sir.  ?s?sl>”*^,:xs"t.7S,;:tSo':ti^|I1?^o.;.3t^..^oo..^ps.. 

Jlfl7sl9?nOo!??l?s!l?2So)^?375.9?SoL.?S?S.,27S0.9?BH.9 

33000. » 312G , »3?B0.,337S. ♦ 3S0  0 • / 

CALL  PLOTS(  0.0,?) 

WPITF (6.7000) 

XSCALF  = 10. /20, 

YSCALE  = 3.  /3000. 

RE *0 (5» 1 03)  NOTS 
WRTTF (6.1 03)  NOTS 
00  8RR  Ka], NOTS 
no  IS  J=1  «6 
no  20  1=1. iso 

LGTl  ( J»  T ) =VFLl  ( J»  I ) = 0. 

0 rOMTTNUF 

S PFAO (S?1 OS)  TITLF.rPATRS 
WR I TF (6.10S)  TITLE* I°AIRS 

PFA^*(S,  110)  TAPAIR.TTTLE.1  *r'*^pJl<'co»rp 

toR  JTE(S.IJI)  lAPAlR.TITLEl.NMPT.SPACF 

WR  T TF  ( f>  • 1 4 0 ) ( PHI  ( I ) * I =1  »NMpT  ) 

MAXWL ( T APA IR)  = M AX toV ( I AP A I R ) = 0. 

MINW(  (IAPAIP)  = MINtoV(IAPAlR)  = 9909, 

no  l oo  T = 1 * nmp t 

**•  CAICULATE  WAVELFNGTH  *** 

TF ( PMI ( T ) . EQ» 0 , ) GO  TO  A01 
WAVLGT ( T ) = 360  * S°ACF  / PHl(T) 

GO  TO 
1 roMTTMJF 


*NOTS  = 
*IAPAIR  = 
•TPAIRS  = 

•titlf  = 
•titlfi  = 
*nmpt  = 
•SPACE  = 
•FRFQ  = 
*PH  I = 
•WAVLGT  = 
•WAVVFL  = 
**#••****< 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 


MiiMPFP  OF  TEST  STATIONS 
number  USFO  TO  identify 

NUMRFP  OF  ACC  PAIRS  AT 


MQORE  MAY 


ON  A RUNWAY 
WHICH  PAIRS 
1 PARTICULAR 
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c 


c 


c 


c 


WAVLGT  ( n = 0. 

40?  CONTINUE 

TP  ( WAVLGT ( T ) .1.  r.MINWI.  ( TAPAIH)  ) m T NWL  ( I Au A T R ) 

TP  (WAVLGT  ( I ) .GT.mAXWI  ( I A° A I R ) ) m AX  «iL  ( 1 ArjA  I R ) 

***  CALCULftTF  W A VP  VFLDCTTY  *** 

WAVVFL(T)  = PWPQ  < I » * WAVLGT(T) 

TF  (wfcVVPL  ( I > .GT.R0O0. > WAWFl.  ( T ) =80  00. 

TP ( WAVVFL U ) .GT.MAXWV ( I APAIH) ) MAXWV ( I APAIR) 

TF  (n'AVVPL  (I  ) .LT.v«INWV(  IAPAIR)  ) MTNWV  ( T A»A  TP) 

LGT1  < IAPATR.  T > = WAVI.GT  ( T ) 

VEL1 (TAOATR.T ) = wAVVFL(T) 

100  TONTTMUP 

***  DRAW  AXPS  FOR  PI  OT  *** 

CALL  AXTS(0.»l..n.»20.,.S.l.  . 2 . . 4HF4 . 0 . 0 . , 1 OhhavFLENGTH,  10) 

CALI  AXTS ( 0.  . 1. , 0.  . 0000. , . OOl. SOO.  . 1000. .4HFS. o.ro. , 1 1HWAVF  vfloct 


WAVLGT < 1 ) 
WAVI.GT  ( 1 ) 


WAWFl  ( T ) 
WAWFL  ( I ) 


1 TY  * 1 3 ) 

CALL  gymroi  (6.O.R.S.O. 14.TTTLP U ) .0.0,20) 

PALL  Symbol (6.0.R.2S.0.14.TITLP (3) .0.0,20) 

PALL  SYmrol(S.(),R.0.0.14,TITLF1  (1)  *0.0.20) 
no  SOO  T = l.NMPT 

WVPLT(T)  = YSCALP  * (WAVVFLID)  ♦ 1. 

TF (WAVLGT ( I) .GT.20. ) W A VLGT ( I ) =20. 

WLPLTII)  = XSCALF  * ( WAVLGT  ( T ) ) 

***  PLOT  WAVE  VFLOCTTY  AS  A FljMPTTOK<  OP  wAVFLFNGTH  *** 

PALL  SYMBOL  ( WLPLT  ( T ) .WVPl.T  ( I ) . . 1 0 . .1 . 0 . 0 . - 1 ) 

SOO  '"ontinijf 

PALL  PLOT  (IS..0..-3) 

***  TARUI  A TF  RESULTS  *** 

WPTTP (6,200) 

WPTTP (6.210) 

WP  I TF  ( 6 » 2 T S ) 

WPITP  (6.220)  ( (I.GT1  ( J,  I ) ,VFL1  ( J.T  ) .J  = 1.6)  , FPF  0 ( I ) .I  = 1«NMPT) 

WPTTP (6,230) 

WPITF  (6.240)  (MAXWL  (J)  .MAXWV(J)  »J  = 1.6) 

WPTTP (6.2S0) 

WPTTP (6,260) (MINWL (J) .MINWV (J) ,J  = 1.6) 

10  PONTIFF 

TF ( IPATRS.PO. 1 ) GO  TO  H«« 

PALL  AXIS(0.,1.»0.»20.«.S»1.»?. .4HF4.0 ,0. , 1 OHW A VFLFNGTH , 1 0 ) 

PALL  AXTS(O..l.,O.,HOOO.,.O01,S0O. . 1 0 0 0 . » 4HF 6 . 0 ♦ 90 . , 1 3HW A VF  VFLOPI 
1 T V • 1 3 ) 

PALL  SYMROL (6.0.S.S.0.14.TITLF (1 ) .0.0.20) 

PALL  SYMBOL (6.0.R.2S.0.14.TITLP (3) .0.0,20) 

PALL  SYMROL (6.O.B.O.O. 14.9HCOMPOSITF.0.0.9) 
no  600  KK=1,TPA[RS 
• IsSAVF  (KK) 
no  610  T=1.NMPT 

WVPLT(I)=YSCALF*(VFL1 (J.T) )♦!. 

WLPLT ( I ) =XSPALF* (LGT1 ( J, I ) ) 

PALL  PLOT  ( WLPI.  T ( T ) .WVPLT  ( I ) . 3) 

PALL  PLOT (WLPLT ( T ) .WVPLT ( I ) .2) 

610  PONTTNUF 
600  ponttmuf 

PALL  PLOT (1S..0..-3) 

BBR  roNTIMIP 

PALL  PLOT (0..0..999) 

STOP  0 

103  FOPMATUM 

IDS  FORMAT (4A1 0.3X. T?) 

]06  FORMAT (//2X.4A10.3X.T?) 

1)0  FORMAT (T2.3X.2A10.TS, PS. 0) 

1 11  FORMAT (?X»T?»3x«?Ain,TS,PS.O) 

140  PORMAT(4F10.0) 

?00  FORMAT  (}H),hx,*1/3*,17x.*2/«.*.)7x»*3/6*.1  7*  .*4/6*.  1 7X , * 1 /S* , 1 7 X , *2 
1/6*) 


e> 
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210  FOMMAT  (IH0.2X«*I  RTfl)  *.4  X ,*VC'I  (1>*.4X.*IGT(2>  (2)  *.4X  • *LGT 

1 (3)  *.4X«*V*l  ( 3)  *.4X  **I_GT  (4)  *.4X  ,*VF|  ( 4 ) * . 4 X , *I_GT  < S ) * . 4 X » *VFL  (*>)  *• 
R4X«*LGT  (ft)  *.*X.*VFL  (M  *,4X»#FPFO*) 

?1S  F O^M AT ( 1 H , ) 

P20  FORMAT n 3F1 0. 3) 

230  FORMAT (1h  , 1 x,*MAX#L  < 1 >*«?X.*MAXWV ( 1) *.2X.*MAXWL (2) *.2X»*MAX#V<2>* 
l.?X,*MAXwl  (3)«.?X,*MAXWV(3>«»?X.*MAX*L(4)*.PX.*MAXrfV/<4>*.?X.*MAXWL 
2(S)*.2X.*MAXrfV(S>*«?X.*MAXiXL(M»»?X.#MAX9V(M*) 

?40  FORMAT ( 12F1 n. 3) 

2-30  FORMAT  ( 1H  . 1X.«MIN*L  ( 1 ) *.?X,*MTMWV(  1 )*«?X.*MIN^l  < 2 ) * « 2X  , *M  J NWV  ( ? > * 
1 ,?X,*mInw|  (3)  *,2X,*MTN*v  (3)  *»?X.«MINWL  ( 4 ) * « 2X  * *mI  mwv  ( 4 ) * . ?X  . *MlNWL 
2(S)  *,2X.*MTMWV  (S) *,?X.*MINWL (M  *.2X ,*MINWV (6) *) 

2*0  FORMAT ( 12F10.3) 

7000  FORMAT ( » 1 * ) 
fno 

SORROliT  T mf  aXTS(X.Y,V1 .VF.SCAI  F , T I C . OLHL  . FMT . THF  T A .LHL  • NC  ) 

C . 

C.  I ATFST  CHANGE  — A (JOUST  7 * 1973. 

C 

OATA  OTW.T^L ANK/0.0174S3?9?b. 1H  / 

cthscos  ( thft  a*oto) 
tcc=o. o**cth 

0TH=S1N ( ThFTA*OT*> ) 
tcs=o.o**sth 

Mr ( VF-VT ) /TIC  ♦ 1.000001 

CALL  PLOT ( X » Y » 3 ) 

M = 0|  UL/TTC  ♦ 0.000001 
X X = X 
YY  = Y 

XINC=CTh*TTC*SCA|  F 
Y JNC=STH*TIC*SCA|  F 

oo  10  1 = 1. M 

F=1 ♦ (m-MOO ( T-l ,M) ) /M 
CALL  PLOT ( X X » Y Y . 2 ) 

FALL  °LOT (XX-F*TCS»YY>F*TCC»?) 

CALL  PLOT(XX*F*TCS»yr-F*TCC«?> 

CALL  PLOT < XX.YY.2) 

XXaXX*XTNC 
10  YY= Y Y ♦ Y 1 NC 

CALL  PLOT  ( X>CH*  ( VF -vi  ) *SC4L.F  • Y^STH*  ( vF-VI ) *SCALF ,2) 

CALL  PLOT ( X « Y • 2) 

XINC=CTM*nLRL*SCA(  F 

yinc=sth*olhl*scalf 
M=  ( VF-VT  ) /OI.9L  ♦1.000001 
OFLXa-STH*SIGN (O.SO.THFTA) -0. 30 
OELY=  CTH*SION(0. 1S,THFTA)-0.0S 
IF  (THETA  .EO.  0.)  DFLY=-.?0 
O0  ?0  1 = 1 .N 
f = T-1 

V = V I ♦F*DLHL 

XX*X^E*X IMC+OFl X ♦ .?* 

YY  = Y*F*YlNC*nF!.Y 

20  FALL  NUMRFP ( XX.YY.O. 1 O.V.O. ,-l ) 

IF  (NC  .LT.  11  * A NO.  LRL  .FO.  ISLAM*)  RF TURN 

XX  = O.F*SCAI_F*  ( VF-VT  ) -0. 0**NC 
yy  = S!F,n<1  ./(?.- AHS(STH)  ) .THETA) 

IF  (ThFTA  ,FO.  0.)  YY=-0.S0 

FPS  = 0.001  f IF  (AHS(THFTA)  .OT.  3*0.)  FPS=-0 . 001 

MS ( 1 AHS( IWT (THFTA-FPS) ) ♦90) /1«0  * SIGN( 1.  .THFTA) 

ANG=THFTA-1 HO.*M 

TF  (MOO (N.2)  ,NF.  0)  XX=XX^0. 1 ?*NC 

XXXsX. XX*CTH-YY*GTH  4 YYY=Y^XX*STH^YY*CTH 

CALL  9YMR0I  (XXX,YYY.0.1S.LPL»ANGiNC) 

PF  TURN 
fnd 
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APPENDIX  D 

PROGRAM  PREDICT  USER'S  MANUAL 

PREDICT  DATA  CARD  SPECIFICATIONS 
(CERF  VERSION) 


A. TITLE  AND  SOLUTION  TYPE  CARD  ( A8,17A4,1X,A2) 
COLUMN  VARIABLE  DESCRIPTION 


1-8 

TITLE  (1 ) 

AIRCRAFT  NAME (LEFT  JUSTIFIED) 

9-76 

TITLE(2) 

NAME  OF  AIRBASE 

RUNWAY  OR  TAXIWAY  NUMBER 

DATE  OR  OTHER  PERTINENT  DATA 

78-79 

LRWIP 

RUNWAY  TYPE 

T = CHANNELIZED  TRAFFIC  AREAS  (LEFT  JUSTIFIED) 
RA  = NON-CHANNELIZED  TRAFFIC  AREAS 

LOAD  FACTORS  AND  PRINT  CONTROLS  (12 ,1X,F3. 0,2(1  X ,11 ) ,5X ,F10. 0,15) 
COLUWJ  VARIABLE  DESCRIPTION 

1-2 

NUMLAY 

NUMBER  OF  MATERIAL  LAYERS  (MAX  = 10) 
(MUST  BE  GREATER  THAN  ZERO  BUT  .LE.  10) 

4-6 

PS  I 

TIRE  PRESSURE  (PSI) 

0 = DEFAULTS  TO  STANDARD  AFCAN  VALUE 
(USER  VALUE  MUST  BE  BETWEEN  70  AND  350  PSI) 

8 

KPPRE 

PRINT  CONTROL  FOR  AFPRE  OUTPUT 
(ALWAYS  USE  1) 

10 

KPPAV 

PRINT  CONTROL  FOR  AFPAV  OUTPUT 
(ALWAYS  USE  L) 

16-25 

WLOAD 

SINGLE  TIRE  LOAD  (POUNDS) 

1 = DEFAULTS  TO  STANDARD  AFCAN  VALUE 
(USER  VALUE  MUST  BE  .GT.  0.) 

26-30 

LCONCOP 

CONCRETE  SPLIT  TENSILE  STRENGTH  IN  PSI 
( WHEN  MATERIAL  TYPE  IS  CONCRETE .LCONCOP  MUST 
BE  GREATER  THAN  ZERO  BUT  .LE.  1500  ) 

LAYER  MATERIAL  PROPERTIES  CARD 
(F4.0.1X.A1 ,11  ,F8.0,1X,F3.2,5X,F1.0,5X,F5.3,2X,F3.0,2X,F3.0) 

( NOTE-  ONE  C-CARD  IS  REQUIRED  FOR  EACH  LAYER,  THAT  IS  NUMBER 
OF  C-CARDS  MUST  EQUAL  NUMLAY  ) 

COLUMN 

VARIABLE 

DESCRIPTION 

1-4 

THICK 

LAYER  THICKNESS  (INCHES) 

( THE  SUM  OF  THICK  FOR  ALL  LAYERS  MUST 
EQUAL  144  INCHES) 

6 

MATID 

MATERIAL  IDENTIFICATION  CODE 
C = PORTLAND  CEMENT  CONCRETE 
A = ASPHALT  WEARING  SURFACE 
B = BASE  OR  SUBBASE 
S - SUBGRADE 
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7 


r 

r 

r 

c 

r 

r 

c 

r 

c 

r 

r 


c 

r 

c 

c 

r 

r 

r 

c 

c 

r 

c 

c 

r 

c 

r 

r 

c 

c 

r 

c 

c 

r 

c 

c 

c 

c 

c 

c 

c 

c 

r 

c 

c 

r 

c 

c 

r 

C 0. 

c 

r 

c 

c 

c 

c 

f 


H — 1 1> 


X = SPECIAL  MATERIAL 

PA  I C HSR  MAItRIAL  PROPERTIES  CHECK 

0 = PRUPTT(l)  MOgT  FALL  *•  I l H 1 N AECAN 

ACCEPTABLE  RANGE 

1 = PERM  I Ij  PROPTY ( 1 | VALuto  uoloIOL 

AFCAN  ACCEPTABLE  RANGE 

( AECAN  ACCEPTABLE  RANGES  E JR  PRuPTY(l)  ) 

ASPHAlI  - .GE.  20000  BUT  .LE.  2000000 
— CONCRETE  - . Gfc . 2000000  OUT  «Lt«  SOOOOOO 
^ BASE  OR  SOHBASE  - .GE.  SOOO  bUT  .LE.  lOOUOu 
SUBGRADE  - .Gt . SOO  BUT  .LE.  20000 
SPECIAL  MATERIAL  - .GE.  SOO 

PROP  T Y ( i ) CORRECTED  PLASTIC  MODULUS  P KOM  THE 
NDPT  VAN 


1 7-10 


31-3S 


38-AO 


PROP  f Y ( 2 ) POISSONS  k A T 1 0 

0 = DEEAULTS  TO  STANDARD  AECAN  VALUt 
(USER  VALUE  MUST  BE  BETWEEN  0.  AND  .b  I 

STYPt  LAYER  SOIL  TYPE 

1 = NON-PLASTIC  WITH  FINES  AND  LOw  PLASWCIfY 

2 = HIGH-PLASTICITY. LL  ABOVE  bO 

3 = CLEAN  SANDS 

A = CLEAN  GRAVELS,  POORLY  GRADED  SAND/GRAVEL 

mixtures 

S = wELL-GRADED  SAND/GRAVEL 

(WHEN  MATtRIAL  TYPt  IS  ASPHALI  OK  CuNLRE • E * 
STYPE  MUbI  EQUAL  0.) 


ERATIO  LAYER  VOID  RA T 1 0, EXPRESSED  AS  A DECIMAL 

(WHEN  MATERIAL  TYPE  IS  ASPHALT  OR  CONCKElE, 
ERAT 10  MUST  EQUAL  0. ) 

(WHEN  MATERIAL  TYPE  IS  h »g»UR  X * ERA  I 1 U 
MUST  RE  BETWEEN  .2  AND  ?.) 

SATUR  LAYER  SATURAT I ON » E XpRESStD  Ao  A PERCENT 

(WHEN  MATtRIAL  TYPE  IS  ASPHALT  OR  CONCRETE, 
SA lUR  MUS I EQUAL  0.  ) 

(when  maIerial  type  io  r»j»uk  x,  jaiur  musi 

BE  BETWEEN  0 AN+  1001 


43-AS  PI  layer  PLASTICITY  index, EXPRcSSEU  as  a PERCENT 

(WHEN  MATERIAL  TYPE  IS  ASPHALT  UR  CUNCktlt.  PI 
MUST  EQUAL  0.) 

(WHEN  MATtRIAL  TYPE  Io  B»S»  OR  X,  PI  MUS I BE 
BETWEEN  0 AND  100) 


end  card  (A3> 

COLUMN  VARIABLE  DESCRIPTION 

1-3  TITlE(I)  COLUMNS  1-3  MUST  HAVE  END  IN  THEM 


(NOTE-ONLY  ONE  U-CARD  IS  RLQUlRtU.  IT  WILL  ALWAYS  bt  THt  lASI 
CARD  IN  THE  DATA  DECK . ) 


100 


n.i.i  mi  mwm  l m 


C 

c 

c 

r 

C 

c 

c 


REPtAT  A1B1  AND  L CARDS  FOR  EACH  PROBLEM 


SAMPLE  PPOHLt  mS 


COLUMN  NUMBER 


1 

5 

in 

20 

30 

40 

Fill 

weak  asi 

phal  r 

-concrete 

4 

1 1 

30000 

4 

A 1 

1 50000 

4 

R 1 

38700 

.30 

5 

.22 

8 

B 1 

32100 

.35 

5 

.25 

128 

R 1 

15000 

.40 

1 

.45 

80 

1 3 

Fill 

STRONG  - 

ASPHAlT-CONCREIE 

4 

1 1 

1 3000 

4 

A 1 

300000 

4 

B1 

88700 

.30 

5 

.22 

8 

B 1 

32100 

.35 

5 

.25 

128 

B 1 

15000 

.40 

1 

.45 

HO 

1 3 

B 52 

asphalt 

TtsT 

PROBLtM — MUL  1 I 

WHEEL 

4 

1 1 

1 

8 

A 1 

721600 

.38 

4 

B 1 

88700 

5 

.22 

8 

B 

32100 

.35 

5 

.25 

124 

B 

13000 

1 

• 45 

80 

13 

50 


60 


70 


80 

KA 


RA 


KA 
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ABBREVIATIONS,  ACRONYMS,  AND  SYMBOLS 


elastic  modulus  (Young's) 
force  vector 

force  at  node  i which  produces  a displacement  in  the  y-direction 
shear  modulus 
Rayleigh  wave  velocity 
shear  wave  velocity 

distance  between  compared  accelerometers 
frequency 

tire  pressure  defined  over  tire  length 

acceleration  constant 

number  of  terms 

half-period  of  function 

phase  velocity 

unit  weight  of  material 

wavel ength 

Poisson's  ratio 

boundary  stresses 

phase  angle 


INITIAL  DISTRIBUTION 


Hq  USAF/PRE  1 
AFISC/PQUAL  1 
Hq  AFSC/DE  1 
Hq  AFSC/DLCAM  1 
Hq  TAC/DEE  ] 
Hq  SAC/DEE  1 
Hq  AFLC/DEE  1 
Hq  ADCOM/DEE  1 
Hq  AUL/LDE  1 
Hq  AAC/DEE  1 
AFIT/Tech  Lib  1 
AFIT/CES  1 
Hq  USAFE/DEE  1 
Hq  PACAF/DEE  1 
USAFA/DFCE  1 
ASD/DEE  1 
SAMSO/DEE  1 
ESD/DEE  1 
AFATL/DLOSL  1 
Hq  ADTC/DEE  1 
Hq  ATC/DE  1 
RADC/DEE  1 
USAF  Reg  Civil  Eng  2 
Chief  of  Eng/DAEN-RDM  1 
AFCEC/EMR  5 
CERL  1 
Hq  MAC/DEMP  1 
Dir,  USAF  Eng  WW  Exp  Stn  1 
AFCEC/PR  1 
AFCEC/EMP  TO 
CEL/Naval  Const  Battalion  Ctr  1 
DDC  12 


107 

(The  reverse  of  this  page  is  blank) 


